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Lecture 1: INTRODUCTION 



The study of integrable models has a long and rich history in condensed 
matter physics, beginning with Bethe's solution of the one dimensional Heisen- 
berg Model [1] and extending to our days when a variety of soluble models 
provide the paradigms that form much of our physical intuition. 

Integrable models are typically defined on a line; they may thus de- 
scribe the physics of actually linear systems such as organic conductors 
or the physics of higher dimensional systems where a particular rotational 
mode has been isolated, as is the case with impurity models. Higher di- 
mensional models, though not integrable, often exhibit properties found in 
lower- dimensional relatives which can be studied non-perturbatively. This 
provides guidance to strongly coupled, often inaccessible, hamiltonians re- 
quired to describe properties of the new materials such as heavy fermions 
systems and High- Temperature Superconductors. 

Rather than offer an overview of the theory of integrability and a list of 
solved models, we have chosen to present here the detailed solutions of two 
models (with an occasional side glance at others) which play an important 
role in modern condensed matter physics. The first is the Hubbard hamil- 
tonian, the prime example of a model incorporating strong correlations, and 
thus playing an important role in the efforts to understand Cuperate Super- 
conductors. The second model is the Kondo hamiltonian, describing dilute 
magnetic alloys, whose properties are basic to the understanding of heavy 
Fermions. The models also provide examples of solutions on the lattice and 
in the continuum respectively. 

In the first two lectures we discuss the steps leading from a given (in- 
tegrable) hamiltonian to a set of algebraic equations that encapsulate the 
physics they contain. In the third lecture we solve the equations governing 
the Kondo model, and in the fourth those governing the Hubbard model. We 
discuss the structure and nature of the elementary excitation and their inter- 
actions, and derive the thermodynamic properties of the models by summing 
over all excitations. 

We intended originally to include a fifth lecture to situate the exact so- 
lution in the framework of the Renormalization Group (RG) approach [2]. 
Time constraints, however, will not allow it. Instead, let me provide a brief 
summary. The RG approach provides a systematic way to explore the low 
energy physics of a hamiltonian by successively integrating out high energy 
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modes and incorporating their effects in tlie resulting tlieory. Tliis way one 
constructs a series of liamiltonians all of which have the same physics in the 
infra-red. This renormalization group "flow" may tend to a fixed point hamil- 
tonian which by construction is scale invariant, and typically much simpler 
than the starting hamiltonian. In many cases it is also conformally invariant 
and then can be completely specified in terms of very few parameters [?]. 
When an exact solution is available these parameters can be computed quite 
easily, without resorting to the much more arduous task of constructing the 
RG flow. Once the fixed point has been identified, correlation functions can 
be written down which describe the asymptotic behavior of the model. The 
calculation of the full correlation function from the exact solution, valid on 
all energy scales, is still a major open problem. 

The Hubbard Model. 

The solution of the model was obtained by Lieb and Wu [4] in the form of 
a set of coupled integral equations applying a method due to Yang [5] . They 
also analyzed the equations for the repulsive interaction at half filling, and 
showed that a charge gap opens no matter how weak the interaction. 

Investigating the spectrum of the model, one finds it contains spinless 
charge excitations (in modern terminology refered to as holons ), as well as 
spin-| excitations carrying no charge (spinons). In the repulsive case away 
from half filling both types of excitations are gapless, while at half filling 
a charge gap opens, inducing a Mott transition at zero coupling. In the 
attractive case a spin gap is always present and the charge excitations are 
gapless. 

In the low energy limit the excitation spectrum becomes linear. This 
fact allows a simple description of the low energy properties of the model: 
the fixed point hamiltonian to which the system flows at long distances is 
expected to be conformally invariant. 

An important property of the model with, perhaps, analogues in higher 
dimensions is the charge-spin decoupling. One finds that charge and spin 
degrees of freedom propagate with different velocities and are characterized 
by different couplings. Thus an electron added to the system will not survive 
as a single particle and hence the fermion correlation function will not possess 
single particle poles (i.e. Z=0), signaling a breakdown of the Fermi-liquid 
assumptions. 

A related, but not identical, phenomenon of complete decoupling occurs 
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in the low energy limit. What one observes is that the charge and spin 
excitations do not interact at all and can be treated separately. (A simple 
way to understand the complete decoupling is to consider a continuum version 
of the model which would describe it at weak coupling and low energies. 
Away from half filling one obtains the g-ology model [6], another integrable 
model, solved under the name of the Chiral Gross-Neveu model [7] [8]. In the 
continuum model the decoupling occurs on all scales, and is due to simple 
kinematical considerations.) This decoupling may be an important clue in 
understanding how a Fermi liquid is destroyed; Already in the free field 
theory there is a complete charge-spin decoupling (however, both degrees of 
freedom are chracterized by the same velocities) in the Bethe-Ansatz basis of 
the Hilbert space, obtained by taking the U ^ limit in the equations below. 
The conventional Fermi liquid description, on the other hand, corresponds 
to the choice of the occupation number basis (Fock basis) in the Hilbert 
space. Both descriptions are valid. When U is turned on, the charge and 
spin excitations still maintain their separate identities, but are modified. 
They acquire different velocities and couplings, they interact with each other 
and among themselves, and the Fermi liquid picture breaks down. Even in 
the low energy limit the spectrum can no longer be described in terms of 
fermionic quasi particles. Put differently, the model flows (in the sense of 
Renormalization Group) to a non Fermi liquid fixed point, the Luttinger 
liquid [9] (a particular c=l conformal field theory.) 

The Hamiltonian is well known 

L L 

H = -^E(C^Hia + h.C.) + Uj2ni^ni^. (1) 

i=l i=l 

The fermionic field tpia annihilates a particle with spin component a at 
site z on a chain with L sites. The first term in the Hamiltonian describes 
hopping from site i to site z + 1 and back, while the second term is a crude 
approximation to a Coulomb repulsion. 

The model has an obvious U{1) X SU{2) symmetry, 

i^ia > e'^'lpia 

i'ia ^ Uabi'ib 

expressing the charge consevation and invariance under spin rotation. The 
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associated generators are given by the spin operators, 
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s. = ^E(«^T - = E^^T^a, s- = {s^r (2) 

^ i=i i=i 

and the number operator, 

L 
i=l 

The conservation law associated with the U(l) symmetry Af allows us to 
study the hamiltonian for a fixed number N of electrons. We shall label the 
states with the quantum numbers M and M' = N — M oi the down-and up 
spins, \F{N — M^M) >, and the corresponding energies E{N — M^M; U). 
The z-component of the total spin is S;, = |(M' - M) = - 2M). By 
construction, this is also the value of the total spin S ^ S = S^^ since the B- 
operators we shall use in Lecture 2 to build the eigenstates commute among 
themselves and thus specify a Young-tableau with definite transformation 
properties. The states we construct are therefore SU{2) highest weight states 
and the rest of the multiplet is obtained by repeated action of the lowering 
operators S~. Beyond N and M, each state is labeled by an infinite set of 
quantum numbers which we shall specify later. 

There is another, less obvious, charge SU(2) invariance (of which the U(l) 
is a subgroup) present in a slightly modified version of the model [10], 

^' = ^-^EKi + «^,T)- (4) 
^ i=l 

We added a chemical potential term to the hamiltonian. In a grand canon- 
ical ensemble the model will be half filled. Equivalently, the symmetry will 
show up if we work in the canonical ensemble and choose the filling appro- 
priately. The symmetry is realized by number density and pair creation and 
annihilation operators, 

C. = \ E(r^iT + ni,) - |, C+ = E(-1)V^t^,,, C" = (C+)*. (5) 

i=l i=l 

As the number operator does not commute with and the symmetry mani- 
fests itself only upon comparing excitations in systems with different number 
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of electrons. Still we shall find some consequences of the symmetry even 
though we mostly work with a fixed number of particles. 

We shall discuss the repulsive as well as the attractive model. The fol- 
lowing particle-hole Z2 transformation 

^;t-(-iMt (6) 
^a-(-im, ^;t-(-iMt (7) 

leads to the relation [4], 

E[N - M, M; U) = [N - M)U + E{N -M,L-M; -U), (8) 

between the the energies of states in the two cases. The eigenstates, we 
shall see, are related in a more complicated way. The eigenvalues of the 
the modified hamiltonian are particle-hole symmetric at half filling and in 
particular one has 

E'{N - M,M-U) = E'{N - M,N - M--U). (9) 

We turn now to study the hamiltonian in a Hilbert space TYjv of N par- 
ticles, defined with respect to the vacuum state |0 > containing none, 

^ia|0 >= 0. (10) 

The states that span TYjv are of the form 

N 

\F>= E i^a,...a^(r^l...r^^^)^Cn.|0>, (H) 

ai...ajvni...njv i=l 

and the Fock eigenvalue problem 

H\F >= E\F > (12) 
turns into its iV-particle version 

hF = EE, (13) 
with the first quantized hamiltonian, 

N 

h = -tJ2Aj + Uj2S^^^,, (14) 
i=i j<l 
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acting on the wave function -Fai..ajv(™i---™iv)- The hopping operator A (the 
discrete version of the Laplacian ) is given by: 

AjFa,..ai^(ni...nj...nN) = Fa,..ai^(ni...nj + I.-.tin) + Fa,..ai^(ni...nj - I.-.tin) 

(15) 

Let us note that when a model is studied perturbatively the starting point 
is the hamiltonian 

Ho = -tJ2ii^:^iai^i^ + h.c.), (16) 

i 

whose ground state is the Fermi- sphere. Subsequently, the interaction is 
turned on and the energy levels are corrected order by order. In some cases, 
when nothing dramatic (such as a phase transition) happens, this procedure 
leads to |0 >, the true ground state. 

Our approach here is different. In the presence of a finite volume (infra- 
red) cut-off X, and in the presence of an ultraviolet cut-off (the lattice spacing, 
in this case) the true ground state and the empty state are in the same Hilbert 
space. Therefore one can use the representation of tp and tp* as creation 
and annihilation operators with respect to |0 > to construct a full set of 
eigenvectors and, in particular, determine |0 > . 

We proceed, then, to diagonalize h within TYat- We shall do it for finite 
iV, and eventually take the thermodynamic limit: N ^ oo,with n = 
N/L fixed. 

Begin by considering the case N=l. Now 

h = -tA (17) 

with the obvious solution: 

Fa{n) = A^e'*'", E = -2t cos k 

Proceed to consider the case N = 2. Then, 

h=-t{A, + A2) + US^,^,. (18) 

The particles interact only when rii = n2 = n. Away from this boundary 
the Hamiltonian is free, and the wave function is given as a product of single 
particle solutions: 

i^ala.(r^l,r^2) = Ae'^^^^+'^-^-[Aa,aAni " ^2) + 5,,a.^(r^2 - r^l)] 
= e*'^^+'*'^'^M^aia.^(rii - n^) + B^,^J{n2 - n^)] 

-e*'^^+''^'^1Aa.al^(r^2 - r^l) + B^^aAni - ^^2)] 
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where A is the antisymmetrizer, and 0{n) is a step function. The correspond- 
ing eigenvalue is 

£; = -2^(cosA;i + cos A;2). (19) 
We introduce now the S'-matrix relating the amplitudes in the two regions 

Baia2 = 'S'aja2^feife2 • (20) 

This is the bare S-matrix, later we shall also discuss the dressed or physical 
S-matrix. To determine it we have to impose two conditions: 

(1) Uniqueness^ to ensure that the value of F{n,n) is defined indepen- 
dently of the region. This leads to 

-^nin2 ^n^ni ^nin^ -^n^ni (^f ) 

in other words 

I-PS = S-P (22) 
where P is the spin exchange operator, and / is the spin identity operator 

{P-^)aia2 — Paia2-^i'lb2 ~ -^a2ai 
(-^^)aia2 — -^ala2^felfe2 ~ -^aia2 

Clearly Pj^lf^ = S^lS^l, and = S^lS^l- Expression (22) can be rearranged 
to 

^(1 + ^)^^(1 + ^) = ^(1 + ^), (23) 

indicating that S' = 1 in the symmetric channel. This is not surprising since 
an on-site interaction operates only in the anti-symmetric spin channel. The 
S-matrix therefore is of the form 

S = ^{l + P) + ^{1-P)s (24) 

where s is a scalar in spin space. 

{2)The Schrodinger equation on the boundary = n2. That is, 

-t[F{n + 1, tc) + F{n - 1, tc) + F{n, n + 1) + F{n, n - I)] + UF{n, n) = EF{n, n). 
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Explicitly 



= -2^(cos ki + cos A;2)e'(*'i+*'^)"^^(/ - P)A. 



Solving for s we find 



where 



z(sin ki — sin k2) + f 
z(sin ki — sin A;2) — f 



(25) 



(26) 

and the S'-matrix S'a^a, becomes (when particles j and / interact), 

oil — objbi _ (sin kj — sin kl)Iaja, + i^Paja, 

~ ~ (sin%-sinA;0 + if ' ^ ' 

Let us generalize the construction to N particles. Begin by dividing the 
configuration space into N\ regions according to the ordering of the parti- 
cles on the line, and label them by elements of the permutation group. For 
example, the region (tcs < rii < n^...) will be labeled by the permutation 
Q = {Qi = 3,(^2 = ^tQs = 5...)eS'jv. As there is no interaction in the 
interior of these regions the wave function will given as a sum over a prod- 
uct of single particle wave functions. In our case F = PLANE WAVES. 
When a boundary is crossed, two particles interact (note that mutliparticle 
interaction is forbidden by Fermi statistics ) and hence the amplitudes in the 
regions across the boundary will be related by the two particle S-matrix just 
determined. We thus consider wave functions of the Bethe-form (The Bethe 
Ansatz): 

Q 

with the energy and momentum given by 

E = ^-2^cos% (29) 

i 

P = Y.kr (30) 
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In eq(28) the Q-sum. runs over all the N\ regions, 0{xq) is equal to 1 if the 
particles are ordered according to Q and vanishes otherwise, and ylai..ajv('5) 
is the spin amplitude in region Q. The amplitude in region Q is related to 
an adjacent amplitude Q\ differing from it by the exchange of neighboring 
particles i and j, via the S-matrix S'*-^, 

-^ai_..ai..aj..aff{Q') = ('S'*'^ )ai'...a^^^l -^iV l*?) = ('S'*'^ )a;aj. ..ajv 5 l*?) 

where for convenience we regard the 2-particle S-matrix S'*-^ as carrying N 
indices operating in an N particle spin space but acting non trivially only on 
particles i and j, 

is'nt:t = n ^ii (31) 

We labeled the S-matrix by the particles it acts on. Let me be more 
explicit: We are considering a region Q = {Ql^ Q2 • • • Ql, Q{1 +!)••• QN) 
and an adjacent region Q' = {Ql, Q2 • • • Q{1 + l),Ql • • • QN). In region Q 
the particle i = Ql is to the left of particle j = Q{1 + 1) while in region Q' 
the particle j is to the left of particle i. In other words, Q' = P^^Q. Hence 
to move from region Q to region Q' we apply S'*-^, A{Q') = S'^-' A{Q) and 
similarly A{Q) = S^'A{Q'). Note that S'^ = (S'-^'')"\ and in the amplitude 
on which S'*-^ acts particle i is to the left of particle j. 

How are regions Qi and Q2 related when they are not adjacent, that is, not 
related by a single transposition? There is always a path (and usually more 
than one) in the permutation group, given as a product of transpositions, 
leading from Qi to Q2- To relate the regions we take the corresponding 
product of the S-matrices. If, for example 

= P^ipi^P^^Q^ (32) 

then 

A{Qr) = S'^S^''S''^A(Q2). (33) 

If the path is not unique, consistency requires that the result be path inde- 
pendent. Consider the N = 3 case. Configuration space is divided into 3!=6 
regions, not all of which are adjacent. We label the regions by (ijk) for the 
ordering (rii < rij < tcj,), and draw a line between adjacent regions. Thus 
one obtains the diagram. 
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(123) ^ > (213) 



(132) 




(312) ^ (321) 



Figure 1: The regions in the Configuration Space of 3 Par- 
ticles. 

Here region (123) is related to region (213) by exchanging particles 1 
and 2. Therefore A(213) = S^^ A(123). Continuing this process one may, 
starting from the wave function in region (123), construct the wavefunction 
in any other region. However, note that region (321) can be reached by two 
different paths and hence for the procedure to be consistent these paths must 
yield the same result. Thus consistency in the three particle case leads to 
the Yang-Baxter Equation (YBE) [5] 

5^235,135,12 ^ 5125,135,23 ^34^ 
in the notation of eq (31), or explicitly 

f C23\d2di f cl3\dib3 j cl2\bib2 _ f cl2\did2 f Cl3\bid3 j c23\b2b3 jr,f-\ 
W Ia2a3\'~' jaidaW Idid2 ~ W iaiaaW JdiajW Jda^ia ' 

where both sides act on 6263 (123) to produce ylaia2a3 (321) . 
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Similar considerations in the case of N=2 and N=A require that the 
matrices also satisfy 

S'^S^' = I (36) 

Furthermore, it can be shown that these relations are sufficient to guarantee 
consistency, that is, path independence, for all N [11]. 
We may therefore rewrite eq(28) as 

Fa,...a^{nr...nN) = Ae'^^''^''^J2S{Q)Aa,...aJ{nQ) (37) 

Q 

with Aai...ajv being the amplitude in some reference region labeled /, Q is 
given as a product of successive transpositions, and S{Q) is the corresponding 
product of S-matrices. 

To Summarize: if the S'-matrix derived from the Hamiltonian satisfies 
the YBE, then the Bethe-form of the wave functions is consistent and the 
model is integrable. A direct calculation shows that indeed the S'-matrix we 
derived for the Hubbard model does satisfy the YBE. 

We may stop at this point to ask how a model can fail to have eigenstates 
of the Bethe-form. The answer is that we insisted that the single particle so- 
lutions whose product forms the wavefunction in the interior of each region 
be labeled by the same set of momenta kj. This goes beyond energy and 
momentum conservation which place the milder requirement that J2j cos kj 
and kj be the same in each region. The fact that the eigenstates can 
be built with momenta kj^ individually preserved under the interaction, is a 
special feature of an integrable model and refiects the fact that the model 
possesses additional dynamical symmetry, expressed by an infinite number of 
commuting conserved charges. A familiar example in quantum mechanics is 
the hydrogen atom where the conserved Lenz vector is responsible for its in- 
tegrability. The existence of these conservation laws restricts the dynamics of 
the model and has many consequences: There is no particle production in on- 
shell collisions, multi-particle (physical) S-matrices factorize into prodcucts 
of two-particle matrices, the Thermodynamics can be consistently expressed 
in terms of these on-shell S-matrices. We shall not pursue this direction 
here but concentrate on the off-shell consequences of the conservation laws 
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which are manifested by the YBE guaranteeing consistent wave function of 
the Bethe-form. 

Periodic Boundary Conditions: To properly quantize the model one needs 
to introduce a volume cut off to regulate the infra-red behaviour. We shall- 
study the model on a finite ring with length L : 

Fa,...aA-nj = 0...) = Fa,...aA-nj = L...). (38) 

This condition, which is imposed in different regions in configuration 
space, can be translated to conditions on the wave function in a single region, 

iZ^t\:XA,,...,JQ) = e-''''^A^,...^JQ), (39) 

where 

To derive this expression consider an amplitude Aai...aj^{Q) and act on it 
first with S-'-''^^ to exchange particle j with particle j ' + 1 to its right, then 
act with S-'-''^^ to exchange it with the particle j + 2 which is adjacent to its 
right after the previous exchange with j + 1, and so on till it is brought to 
the right end of the system. Now repeat the operation, moving it to the left; 
act on Aai...aj^{Q) with S-'~^-'^ then with S-'~^-' and so on, till the particle j 
reaches the left end of the system. From periodicity 

S^^....S^-^^S^-^^A(Q) = S^^....S"+^S"+^A(Q)e''''^ (41) 

leading to eq(39). 

We need to diagonalize the spin operator Zj acting in the N spin space, 
(in other words, Zj acts on spin functions ^6^ . 6^^ G = Hj Vj^ where Vj = 

is the spin space of particle j). From the eigenvalues of the new spin hamil- 
tonian Zj we can find the energy and momentum of the Fock eigenstate \F > 
constructed from the spin eigenfunction Ah^,,,bj^- Note that the Zj commute, 
[Zj, Zi] = as a result of YBE , and can be simultaneously diagonalized. 

The diagonalization of the (spin-Hamiltonian) Zj was achieved by Yang 
[5] by means of another Bethe Ansatz built out of single "particle" spin 
functions. The role of the vacuum now is played by the ferromagnetic state 
(all spins aligned), a single particle state corresponds to a single spin flip and 



13 



a general eigenstate is built out of products of single particle solutions. This 
approach is reviewed in [12] [13]. A related problem arose in the study of the 
6-vertex model where the spin operator Z plays the role of the transfer matrix 
[14]. Algebraic methods for its diagonalization were developed by Baxter [15], 
and further extended by the St Petersburg school under the name of Inverse 
Scattering Method [16]. We shall discuss the diagonalization of Z using the 
latter approach. But before doing so let us turn to the Kondo model and 
construct its eigenstates. We shall find that similar spin problem Z arises. 

The Kondo Model. 

The Kondo model describes the interaction of a conduction band with a 
localized spin impurity. The litterature on the subject is immense [17]; here 
we concentrate on some of the theoretical aspects . 

The conduction band is described by the Hamiltonian, 

iJo = E^(^)iA (42) 

k 

with c^^ annihilating an electron with momentum h and spin component a. 

The conduction band is coupled via spin exchange interaction to a spin 
(To localized at r = 0, 

El = J*:(f = 0)^,fe*fe(f =0)-ao (43) 

The field ^a{^) is the Fourier transform of c^^. 

Since the model we consider is rotationally but not translationally in- 
variant, an appropriate basis for the electron annihilation operators is Ckim,a^ 
expanded in angular modes around the impurity, rather than c^^. Of these 
modes we assume now that only the s-wave modes have non-zero coupling 
to the impurity. Later we shall discuss the case where higher orbital modes 
couple to an impurity leading to the Multi- Channel Kondo model. 

We further restrict our attention to low energy phenomena, entitling us 
to retain only momenta k close to the Fermi surface: k = kp -\- \q\ <^ D, 
where D is a cut-off, of the order of kp^ which will be considered large 
compared to any other physical scale in the problem. Linearizing the energy, 
e{k) = e^kp) + vpq^ and Fourier transforming with respect to q we find that 
the ^-component of Hq becomes in the limit D — )• oo, 

Ho = -iJ raix)d^M^)dx. (44) 
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We have chosen our units so that vp = 1. The field tpa{x) is the Fourier 
transform of Cfc^+g,oo,a, where x is the variable conjugate to g, — oo < a; < oo. 

We obtain a l-d field theory with only right moving electrons, a reflection 
of the fact that the Fermi surface is simply connected. One may also carry 
out an analogous derivation in real space, where the problem is defined on 
the half line with incoming and outgoing waves, then mapping the outgoing 
to incoming waves but defined for a; < 0, one obtains a full line problem with 
only one kind of movers. 

Adding the interaction term to Hq we obtain the Kondo model which is 
the starting point of our investigation: 

H = -i j raix)d:,M^)dx + JC(0)^afe^fe(0) • ao + JV:(0)^a(0). (45) 

We also included a term that couples to the electron charge density at a; = 
with strength J'. Since we have linearized the model charge and spin degrees 
of freedom completely decouple on all scales in Hq. This is not modified by 
the interaction terms which couple separately to the spin and to the charge 
of the conduction electrons. 

The linearization procedure is valid only when all energy scales (such as 
temperature T, magnetic field excitation energy e) are small compared 
to the cutoff. Otherwise the linearization breaks down, and details of the 
band structure (reflected in the cutoff procedure) become relevant. We shall 
consider only quantities that characterize low-energy properties of the model, 
and are independent of the cutoff scheme. These quantities we shall call 
universal. 

As long as we consider only universal quantities we do not have to insist 
on a particular cutoff scheme. Different schemes may be employed to give 
the same universal quantities, though outside their domain of applicability 
universality may break down and results may vary; if we want to analyze 
properties at T ~ D, much more care must be taken in the construction of 
the model, and there may be only one physically acceptable cutoff scheme. 

The above observations allow us to apply the considerations and methods 
of quantum field theory to the problem. Thus, as the coupling constant 
J is dimensionless, the Hamiltonian is renormalizable, and divergences are 
expected in the calculations. The divergences are absorbed within the cut- 
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off scheme chosen[18], and any particular numerical value assigned to J is 
defined only with respect to that scheme. 

One might think that no scale remains in the problem, as the coupling 
constant is dimensionless and the cutoff is considered infinite if we restrict 
ourselves to the low energy regime. 

One of the fundamental properties of the model, however, is the appear- 
ance of dynamically generated scale To (to be defined later) which uniquely 
determines the low-energy physics. This scale depends on the cutoff D and 
the coupling constant J in the following generic way: To = -D exp[— a/A( J)], 
where A(J)— )-JasJ— )-0. The explicit form of A (J) depends on the scheme 
used. In the conventional momentum cutoff scheme, p < Dm, one finds. 

To = i)^e-("/-^^)+(^/')i"-^^+-. (46) 

This scheme, however, spoils intergrability which is restored only when Dm is 
taken to infinity. We impose therefore another cut-off respecting integrability 
while still finite [19]. In this scheme one finds 

To = De-""'^. (47) 

Still, both constructions are characterized by the choice of, say. To = 
0.0007 eV. This value is the only relevant scale in the scaling (universal) 
regime which defines the low-temperature and low magnetic field properties 
of the model. In this region the free energy F takes the form 



^{T,h;D,J)^f 



T_ h 

Yo'Ti 



h,T<.D (48) 



where the function / is universal in the sense that it is independent of the 
particular scheme used to define the model. The cutoff and coupling constant 
enter only in the combination determining To- Also, any other scale must be 
related to To by pure numbers that are directly calculable. These numbers 
are universal. 

The part of the scaling region where T ^ To will be called the high- 
temperature region (still T ^ D). As we shall see, this is the weak coupling 
regime, where the effective coupling constant is small and the physics can be 
captured by perturbing around the weak coupling fixed point hamiltonian, 
Hq. The low-temperature region (T <^ To), however, is a strong coupling 
regime governed another fixed point hamiltonian, H* describing a local Fermi 
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liquid. The crossover in behavior from the strong coupling regime to the weak 
coupling, can be described as a renormalization group flow in the space of 
effective hamiltonians, and is the essence of the Kondo problem [20] [21] [22]. 

The crossover can be driven by any physical parameter. Let us discuss 
it as a function of the temperature. Consider the impurity susceptibility %*, 
which is the term in the susceptibility left over after subtracting from the 
total susceptibility % = dM/dH the contribution of the electrons. (We take 
electrons and impurity to have the same g factor.) 

As we shall see, the high-temperature region lies in the weak coupling 
regime allowing us to apply perturbation theory to find that the impurity 
susceptibility attains its free value = fJ,^ /T (Curie law) up to corrections 
that vanish logarithmically at high temperatures: 

where a new scale Tk has been defined by requiring that the l/[ln^(r/rfc)] 
term be absent. This is a normalization condition on Tj,, the high-temperature 
or perturbative scale, which is conventionally referred to as the Kondo tem- 
perature. 

While the high-temperature region is thus accessible by perturbation the- 
ory, the system enters a strong coupling regime at low temperatures and its 
properties change drastically. Due to the strong coupling to the electrons the 
impurity spin will be screened, leading to a finite susceptibility, %q, at zero 
temperature. Thus define the scale To, 

x' = ^ 

Trio 

for the low-energy regime of the model. The ratio 

W = ^ (51) 

o 

is a universal number characterizing the temperature crossover. The crossover 
occurs, of course, as a function of other quantities, be it the excitation en- 
ergy, the magnetic field or any other energy scale, each crossover having its 
own universal number. These numbers relate different asymptotic regimes 
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(50) 



and as such cannot be calculated from an effective hamilonian probing only 
the neighborhood of one fixed point. Instead, a complete construction is re- 
quired, valid over all scales. This was first carried out numerically by Wilson 
[21]. 

We mentioned already that the physics of the strong coupling regime 
of the model is Fermi liquid like; it is characterized by a local potential 
center (the remnant of the impurity that was screened), and by induced 
weak interactions among the electrons. We shall discuss now a genralization, 
the multichannel Kondo model, where other, non-Fermi liquid fixed points 
may be reached in the infra-red. 

The model was introduced by Nozieres and Blandin to describe "real 
metals" [23]. Taking account of the the orbital structure of the impurity 
they derived the most general exchange hamiltonian to describe the Kondo 
effect. When the atomic shell (with orbital quantum number /) is half filled, 
Hund's rule indicates that the ground state is an orbital singlet with total 
spin S = (2/ + l)/2. The electrons scattering off the impurity then also carry 
the orbital quantum number m, — / < m < /, and one ends up with a 
multichannel version of the Kondo model, 

H = -i J i>*^^{x)d^i>a,mix)dx + Ji>*^^{0)aabi>b,miO) ■ S (52) 

Here m = l,...,/ = 2/ + lis the orbital channel (or flavor) index and the spin 
operator S is in spin-S' representation of SU(2). In the hamiltonian the values 
of / and S are unrestricted, though in a magnetic impurity hamiltonian 
/ = 25'. Other non-magnetic applications of the model exist with other 
values of spin and flavor [24] . 

The nature of the infra-red fixed point depends on those values [23] : for 
/ < 25' the coupling J flows to infinity leading to a screened impurity in the 
case / = 25', and to a partially screened impurity 5'' = 5' — | in the case 
/ < 25'. The strong coupling fixed point becomes unstable when / > 25' and 
the infra-red physics is controled by a new, finite coupling fixed point. This 
new fixed point is expected to describe non Fermi-liquid behavior. 

We shall show next that the Kondo hamiltonian is integrable [25] [26], 
and construct a complete set of eigenstates. We shall find that the spec- 
trum consists of spin-1/2 uncharged excitations, spinous, as well as spinless 
particles carrying the charge degrees of freedom, holons [25]. The spectrum 
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bears similarities to the spectrum of the Hubbard model we discussed ear- 
lier. In fact, in all integrable models belonging to this class, in a sense to 
be discussed, the fundamental excitations carry the same quantum numbers. 
They may differ in their dynamics, though; unlike the spinous and holons 
in the Hubbard model, these excitations decouple on all scales as result of 
our choice of a linear spectrum. The spinous however interact with each 
other and we shall calculate the physical S-matrix. We shall also derive an 
expression for the phase shift a spinon undergoes as it passes the impurity. 
Then we shall proceed to calculate the free energy, and explore the crossover 
as a function of thermodynamic parameters such as temperature and mag- 
netic field, or dynamic parameters as excitation energy and momentum. In 
particular, we shall find an analytic expression for as well as for other 
crossover numbers. 

We turn now to the diagonalization of the model. In the Hilbert space 
Ti-N^ of electrons the hamiltonian takes the form 

h = -iJ2d.,+Jj2 ^(^i)^i • ''"o + J'J2 ^(^i) (53) 
i=i i i 

A single electron interacting with the impurity will be described by a 
wave function 

Fa,ao(^i) = e"^^nAa,aA-^j) + Ba.aA^j)] (54) 

here aj and ao are the spin indices of the elecrton and the impurity respec- 
tively. This function obviously satisfies the Hamiltonian for Xj ^ and has 
the eigenvalue: 

E = k (55) 
Applying the Hamiltonian and evaluating it at Xj = we have 

(h - E)F(O) = -i(B - A)S(O) + (Jaj ■ ao + J')(1/2)(A + B)S(O) = (56) 

where we chose the convention (or rather the renormalization prescription) 
6{x)S{x) = {l/2)S{x) = 6{—x)S{x). Hence the S-matrix that relates the 
amplitude Bajao to the amplitude Aajao is, 

•0 ^ i + (1/2) J^,- • ^0 + (1/2) J- ^ i + JPi° + J" 

i-(l/2)Jaj-ao-(l/2)J' i-JP^^-J" ^ ' 
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where we used = (l/2)(/ + ■ (To), and set J" = (1/2)(J' - J). The 
S-matrix can be brought to the expHcit form, 

p'o _ ,>pio 

S'^'o = e-^^[a,-„P° + &,-oP^'°] = e-^^— -— — (58) 



with 



2J 



1 + J"2 - J2 
i, ^ 1 + J--^ -J' + 2iJ 

^ 1 - J"2 + J2 - 2iJ" 

When we wish to proceed and construct the iV-electron wave function 
a problem arises: the hamiltonian does not contain any interaction terms 
among electrons. That might induce us to adopt S'*-^ = / as the scattering 
matrix of electron i and j, but this choice would not satisfy the YBE 

In fact the non commutativity of 5'-^° and 5'*° captures some important aspects 
of the model: after electron i crossed the impurity the latter is left ia a 
different state then before. Hence the state in which electron j finds the 
impurity depends on whether it crosses the impurity before or after electron i. 
Herein lies the differnce between a system of electrons interacting with a fixed 
potential (a one-body problem since all electron see the same potential) and 
a Kondo system, where the impurity correlates the motion of all electrons. 

This non- commutativity, however, does not ruin the integrability of the 
model. Considering more carefully the model for two electrons away from 
the impurity 

h = -i{di + dj) (60) 

we note that an arbitrary electron-electron S-matrix may be introduced, and 
one is allowed to consider a basis of free eigenstates of the form 

Fa.a,(xiXj) = Ae'^'"^'+'''^'^[Aa,aA^i " Xj) + (SA)a,aA^j " Xi)]. (61) 

This function obviously satisfies the free hamiltonian with eigenvalue E = 
ki + kj for any choice of S', but is not an eigenstate of the individual momen- 
tum operators unless S = I . Let me elaborate this simple point. One would 
tend to write a solution for in the form 

Fa.aAxiXj) = ^e^(*'-+*'^-^)A,... (62) 
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which is indeed an eigenstate of individual momentum operators. As this 
state is degenerate with 

Fia.i^i^i) = ^e^((^-+^)--+(^^-^)-^)A,.,^. (63) 

for any one may sum over p with appropriate coefficients to form eq (61). 
What one is doing, in fact, is to find the zero order approximation in a degen- 
erate perturbation theory. Physically, the eigenfunction can be constructed 
with arbitrary S-matrix because the particles move with the same velocity 
and never cross. 

We conclude, then, that we are at liberty to choose as a basis of states 
one which is determined by an S'*-^ satisfying the YBE. It is clear that the 
choice 

S'' = P'' (64) 
does so, and thus leads to a consistent solution, with 

i^a,...a^e,ao(^l...^Are) = ^^^^^ ^ < ...,^^,,/(xq ) (65) 

Q 

corresponding to the eigenvalue 

E = J2kj. (66) 
i=i 

In eq (65) Q G Sn^+i describes the ordering of the electrons and of 
the impurity, localized at a; = 0. The antisymmetrizer acts on the electron 
variables only. 

We have thus constructed a consistent Bethe-Ansatz for the model, with 
the S-matrix given by 

{^-i4> 11° -xcPi^ electron-impurity 
(67) 
P-^' electron-electron 

indicating that it is integrable. 

Again, imposing periodic boundary conditions one is led to to the problem 
of diagonalizing the operator Z defined in eq(40) but now constructed with 
the (bare) S-matrices corresponding to the Kondo-model. It takes the form 

/ Tjo_-pjo 
(^i)aii = P"~^...P'^P'^ ...e^^- ^^...P"+^ (68) 



IC j 



■a-N 
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We denoted N = + 1, the number of spins in the problem. 

Unhke the case of the Hubbard model, there is only one Z Hamiltonian 
to diagonalize, so that in a given state all electron momenta will be shifted 
by the same amount. In other words, the phase shift of the electrons due to 
their interactions is independent of their motion. This circumstance is due 
to the fact that the coupling constant J is dimensionless, hence the S-matrix 
cannot depend on the momenta. In the Hubbard model, U is dimensionful 
and the S-matrix is of the form S = S{U/{kj — ki)) leading to a coupling of 
all the modes. In the next lecture we turn to the problem of diagonalization 
of the Z-hamiltonian. 
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Lecture 2: The Quantum Inverse Method 



We wish to find the eigenvalues (and eigenunctions) of the operator 
where 

. ^ (a,-a.)ri + icP^i 
ttj — aj + ic 

This operator occurs in the solution of the Hubbard model where aj = 
sin kj and of the Kondo model where aj = 1 or depending on whether j 
refers to an electron or to the impurity, respectively. The solution of other 
integrable models possessing SU(2) symmetry leads again to the same oper- 
ator, differing only as to the values taken by the variables aj. For example 
in the Heisenberg model aj = [1], and in the Backscattering model Model 
ttj = ±1 depending whether j refers to a left mover or a right mover [2]. 

The natural question to ask is whether the Hamiltonians Zj are integrable. 
That this is the case was shown by Yang [4] who solved it by means of 
another Bethe Ansatz. Let us follow a different path, and ask whether the 
hamiltonians Zj possess an infinite set of conserved charges, whose presence 
would guarantee integrability. 

We have already seen that [Zj, Zj]=0. It will prove very useful to obtain 
a continuous version of this statement, since then we may expand in the con- 
tinuous parameter (the one continuing the discrete index j) and obtain a set 
of charges commuting with Zi. One is led thus, following Baxter [5], to intro- 
duce a continuous parameter, usually refered to as the spectral parameter, 
into the definition of the S-matrix 

, , od I %cf^ , , _ , , , , 

S{a) = —^ = a(a)I + b(a)P, 3 

a + zc 

in a way that a continuous version of YBE is also satisfied 

S'^^ia - ^)S'''ia)S^'iP) = S^'iP)S'''ia)S''^a - (4) 

We proceed now to introduce the monodromy matrix [6]. Its construction 
is natural when starting from the 6-vertex model, although its introduction 
here may seem ad-hoc. To begin with, introduce an auxiliary spin space Va 
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which will help us keep track of the proliferating spin indices, and define an 
S-matrix acting in Vj X Va, where Vj is the spin space of particle j. 

a -\- tc 

the variables m, v live in the auxiliary spin space, and a, b live in the physical 
spin space. Now, the monodromy matrix H(a) is defined by: 

E(a) = S^'^ia - ai)S''^(a - as) 5''^^(« - «Ar) (6) 

Where the aj are the physical values appropriate to the model, and the 
product is carried out only in the auxiliary space. Explicitly 

j'Z:\bi...b]^,v / cil^\fei,si / ci2^\fe2,«2 / cN'A\b]^,v /rj\ 

V'~'/ai ...ajv," / J \ /ai,u W /a2,si W /ajv,«jv-l ^ ' 

S1...SJV-1 

where we suppressed the spectral parameters. It is convenient to represent 
the monodromy matrix graphically as: 



2 *i •••^N'" 
bi . . .bN,v 



bi b2 



bN 





Sl 


S2 


S3 SN-l 















ai 



Figure 2: The Monodromy Matrix. 



where the j vertex represents the j S-matrix with the horizontal lines 
carrying the auxiliary variables and the vertical lines at the j-site earring the 
variables of V-' . 

Define also the transfer matrix Z{a) = ^r^S(a), by taking a trace over 
the auxiliary variables. Then using: S'-^^(o) = P-^^, it follows [3], 

Z{a = aj) = Zj. (8) 
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We can now reformulate our question about the integrability of the spin 
problem. If we can show that for /3 7^ a, [Z{a)^Z{(3)] = 0, then expanding 
in powers of (3 we produce an infinite set of charges commuting with the 
hamiltonian Z{a) guaranteeing integrability. 

A sufficient condition to assure the commutativity of the operators Z{a) 
and Z{(3) is the existence of a matrix R acting in Va X Va satisfying 

R (S(a) X Ei(3)) = x S(a)) R (9) 

namely, 

S(a): E{(3)t = E{(3r: E{at Rlf, (10) 

since taking the trace over the auxiliary variables leads to the desired ex- 
pression. In eq (10) the physical indices are contracted in the usual quantum 
mechanical way. Thus rewriting it in its full index glory we have 

To show that eq(ll) holds it is sufficient to examine it locally since the 
monodromy-matrix is given as an ordered product of the S-matrices. In other 
words, consider the problem where S is built of one vertex only, 

RSxS = SxSR (12) 

namely 

from which eq(9) can be deduced by applying eq(12) repeatedly. Comparing 
eq(13) with eq(4) we see that it is again the ubiquitous YBE if we identify 

R'ul. = S{a-(3)l;%. (14) 

In other words, we have shown that a martix i?, implementing eq(9), does 
exist and is given by (all action is in the auxiliary spaces Va X Va) 

R = S(a-/3) P = — . 15 

[a — p) -\- ic 

As a consequence, we have shown the commutativity of the Z-matrices and 
the integrability of the spin problem. 
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Let us consider the construction more explicitly in the auxiliary space. 
The 4x4 matrix R takes the form 



R{ 



a 



fa + b 

b a 

a b 

V 



\ 





a + bj 



(\ 





a.-\-%c a.-\-%c 



a.-\-%c 

VO 



0\ 





(16) 



a.-\-%c 

1/ 



and the monodromy matrix 



fel...fejV 



(AY 




(17) 



with the operators A, 5, Cand D acting in the physical space . Carrying 
out the products as indicated the fundamental relation eq(lO), and equating 
term by term left and right hand sides we obtain the following algebraic 
relations (we display only those necessary to our purposes): 



D{a)B{(3) 



u{P - a)B{P)A{a) + v{P - a)B{a)A{P) 
u{a - (3)B{(3)D{a) + v{a - (3)B{a)D{(3) 



(18) 
(19) 
(20) 



where 



Furthermore, 



v{a) 



B(a)B(P) 
A(a)A(P) 
D{a)D{(3) 



1 


a -\- ic 


a{a) 


a 


b(a) 


ic 


a{a 


a 



B{(3)B{a) 
A(P)A(a) 
D{(3)D{a) 



(21) 
(22) 
(23) 



We shall use these relations to diagonalize the transfer matrix Z{a), explicitly 
given by 
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^fel...fejV 
ai ...ajv 



(«) = <-:.t(«) + ^a;:::t(«)- 



(24) 



The key point in the diagonahzation procedure is the observation that 
the operator -B(/3) plays a role of a creation operator with respect to the 
Hamiltonian Ai^a) + -D(/3), up to "unwanted terms" (generated by the second 
term on the right hand side of eqs(18,19) , when acting on a ferromagnetic 
"up" vacuum |a; >, 



We may show that this ferromagnetic vacuum is an eigenstate of the 
Z hamiltonian using the fact that the latter is given as a product of local 
vertices S^^ which, when written out in auxiliary space, take the form 



We used the identity {P^^)% = (l/2)(^^^^ + ((Tj)l ■ {(TaYJ to separate the 



physical and auxiliary spaces. Ij denotes the two dimensional unit matrix 




(25) 



5^'^ (a) = {a + b/2){a)ljlA + b/2{a)aj ■ a a 




(26) 



acting in Vj. Acting on 



1 





) 



the vertex S-''^ becomes triangular 




(27) 



so that 




(28) 



and we see that \<jO > is an eigenstate of of A{a) and of D{a) 




(29) 
(30) 



28 



where 

^(«) = n — f-- 

j = l 3 ' 

To proceed and find all eigenstates of Z{a) we consider states formed by 
repeated application of the flipping operator B on the ferromagnetic eigensate 

\Ai(3,...(3M)>=Bi(3,)...Bi(3M)\u;>= E ^h-^M'^J.-'^jJ^ > (32) 

ii ■■■Jm 

where the usual spin amplitude notation Aai...aj^ is written as Aj^^^^jj^ by 
specifying the position of the M down spins. Acting on the state with the 
Hamiltonian, and applying the algebraic relations eqs(18-23) we find that in 
addition to terms of the form of the original state, we also obtain unwanted 
terms preventing the state from being an eigenstate. However these unwanted 
terms can be removed by a proper choice ot the parameters (3i....(3m- To 
illustrate the procedure consider M=2. Then, moving A{a) past B{(3i) and 
B{j32) we have 



{A{a) + D{a))B{l3^)B{l32)\uj > = 

u(/3^ - a)u(/32 - a)B(/3^)B(/32)A(a)\u; > 
+u{a - Pr)u{a - P2)B{Pr)B{P2)D{a)\u, > 

+ [u{(3^ - a)v{(32 -a) + v{(3^ - a)v{(32 - (3^)]B{a)B{(3^)A{(32)\u; > 
+ [uia - f3^)via - /^s) + via - f3^)vip^ - /32)]5(a)5(/3i)D(/32)|a; > 
+v{Pr - a)u{P2 - Pi)B{a)B{P2)A{Pr)\u, > 
+v(a - /3^)u(/3^ - /32)B(a)B(/32)D(/3^)\u; > 
= X{a,(3^(32)B{(3^)B{(32)\co > +X^{a, (3^(32)B{a)B{(32)\co > +X2{a, (3^(32)B{a)B{(3^)\co > 

where 

A(a,/3i/32) = u{(3i — a)u{(32 — a) -\- A{a)u{a — (3i)u{a — (32) (33) 

Ai(a,/3i/32) = v(/3^-a)[u(/32-/3^)-u(/3^-/32)A(/3^)] (34) 

X2(a,/3^/32) = v(/32-a)[u(/3^-/32)-u(/32-/3^)A(/32)] (35) 

The condition for B{(3i)B{(32)\co > to be an eigenstate is that (3i and (32 be 
chosen so that 

X^{a,(3^(32) = 7 = 1,2- (36) 
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The generalization to arbitrary M is straightforward [6]: B{(3i)....B{(3m)\<^ > 
is an eigenstate of Z{a) = A{a) + D{a)^ with eigenvalue 



M M 

X{a,(3i...(3M) = n ™(/^7 - a) + A(a) JJ - /^7) (37) 

7=1 7=1 

if the paramaters (3i....(3m are chosen so as to eliminate the "unwanted terms", 
namely: 

X^{a,(3^....(3M) = 7 = 1....M, (38) 

where. 



X^(a,/3i.../3M) =v(/3^- a) 



M M 



5=1,55^7 5=1,57^7 



(39) 



Recall now that our original goal was to find the eigenalues of the oerator 
Zj = Z{a = aj). Denoting these by Zj, we have (see eq(37)) 

= Ma,,f3,....f3M) = n ^I-^l±l£ (40) 
7=1 Pf '^j 

with the parameters (3i...(3m satisfying 

■j-j /3g — /3-y -|- ic -j-j /3-y — aj (41) 

5=1,5^^7 l^s ~ /3j ~ ic /3j — OLi -\- ic 

We may cast these equations in a more appealing form by changing variables: 
(3^ = — ic/2, and recalling that periodic boundary conditiond impose 
Zj = e"**'-'^, we finally obtain: 



ikjL TT f •? 



^ - aj - ic/2 
- aj + ic/2 

and 

Ma a , • iV 



n ? 7 , : (42) 



„ As - + ic _ „ A^ - aj - ic/2 
sj}^^ As-A^-ic~ l}^ A^-ai + ic/2 ' 
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We have now completed the diagonalization of the Z-operator, and have 
solved at the same time the underlying spin problem for the whole class of 
integrable SU(2)-invariant models built with the i?-matrix in eq(16). 

For the Hubbard model aj = sin kj and the Bethe-Ansatz equations take 
the form [7]: 

., , ^ A-v — sin — ic/2 

^J-^ A^ - sm% + zc/2 ^ ' 



anc 



A5 - A^ + zc -i^ A^ - sin A;^ - zc/2 

n A„ A .- =11 A .^^U. , J/0 (45) 



A5 — A-y — ic f-J^ A-y — sin ki + ic/2 



withc=f = |. 



For the Kondo model aj = 1,0 and the equations become [8] [9]: 



and 



il A, - A-y + ic ^ / A-y -1 -ic/2 A"^^ / A-y- ic/2 \ 
^Jj^-y A5 - A-y - ic l,A-y-l + ic/2j l,A-y + ic/2j ^ ' 



with the coupling c given by 



2J 



1 + J"2 - J2 



(48) 



We note that the equations have decoupled; the spin variables A-y are deter- 
mined independently of the momenta kj^ reflecting the charge-spin decou- 
pling discussed earlier. 

We proceed now to solve the equations, and discuss the physics of each 
model. 
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Lecture 3: The Kondo Model 



We now turn to derive the physics of the Kondo model from the Bethe 
Ansatz equations. We shall mainly discuss the excitation spectrum and the 
thermodynamics, with only a casual discussion of transport properties. 

We show here that the fundamental spin excitations are interacting spin- 
1/2 particles, spinons, carrying no charge; while the charge excitations, the 
holons, carry no spin, and are non-interacting [1]. More complex excitations 
are superpositions of the fundamental ones. The spin and charge sectors in 
the model decouple completely. A physical electron, however, is a coherent 
superposition of states in both sectors. Hence the electron-electron correla- 
tion function has no single particle poles. 

We then proceed to obtain the free energy of the model in terms of coupled 
integral equations [2], [3]. The equations are studied in the utraviolet (high 
temperature, or large magnetic fields) as well as in the infra-red regime, and 
a remarkable change in the impurity properties is observed. In the ultraviolet 
limit the impurity behaves as an almost free-spin, given by an effective weakly 
coupled theory, while in the infra-red it is completely screened, becoming 
a non-magnetic scattering center. The properties of the impurity in this 
regime are given by another effective theory, strongly coupled, describing a 
local Fermi liquid. This crossover from one regime to another, the Kondo 
Effect, can be driven by any energy parameter - be it the temperature, the 
magnetic field or an excitation energy. It represents, in the language of 
Renormalization Group, a flow from from one fixed fixed point to another. 
Understanding the crossover requires the construction of the model over all 
energy scales, which is not feasible by methods valid only in the vicinity 
of one fixed point, such as perturbation theory, strong coupling expansion 
or conformal field theory. The crossover was first carried out numerically 
by explicitly constructing the flow in the space of effective hamiltonians [5] . 
Here we study it in the framework of the Bethe ansatz [6], [7]. 

We shall end this lecture with a brief discussion of the Multichannel 
Kondo model. Again a crossover takes place from a weak coupling regime 
in the ulraviolet to a new regime in the infrared. The latter is typically non 
Fermi liquid. 

Let us now determine the energy eigenvalues. Consider a system of N'^ 
electrons on a ring of length L interacting with an impurity localized at 
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X = 0. There are therefore N = + 1 spins in the problem, with N to 
be taken even. An eigenstate will be also labeled by the quantum number 
M, the conserved number of down spins, as well as by an infinite set of local 
quantum numbers {rij^I^} defined below. The z component of the spin of 
the state is given by = ^(N — 2M). This is also the total spin of the 
state since by construction the B{A^) operators (from Lecture 2) commute 
among themselves and thus specify a unique Young Tableau. The states we 
construct are therefore highest weight states, |S', = S >^ and the rest of the 
the states in the multiplet are obtained by the action of lowering operators 

S-. 

In previous lectures we found that the energy eigenvalues are given by 

^ = E% (1) 

i=i 

with the momenta kj obtained from the eigenvalue z 

., , A 1 + ic/2 

ii - 1 - ic/2 ^ ' 

Hence, 

27r M 

% = — r^,• + -^[0(2A,-2)-7r], (3) 

7=1 

with rij an integer arising from taking the logarithm, and 

©(x) = -2tan-^(a;/c). (4) 
The expression for the energy becomes (dropping inessential terms) 

E = j:-n, + Dj:[@i2A,-2)-^], (5) 
i=i ^ 7=1 

where D = / L is the electron density. 

The spin momenta Ai....Am are found from the condition guaranteeing 
the cancellation of the "unwanted terms" (see Lecture 2), 

^ As-A-y + ic / A^-l-ic/2 Y^ / A^-ic/2 \ 
l}^As-A^-ic~\A^-l + ic/2) \A^ + ic/2)' ^' 
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Note that the equations determining the {A-y} have decoupled from those de- 
termining the momenta {%}, reflecting the complete charge-spin decoupling 
discussed earlier. 

Upon taking the logarithm of eqs(6) we find that the variables {A-y} satisfy 
the following set of coupled equations. 

M 

iV^0(2A,y - 2) + 0(2A,y) = -27r/,y + ^0(A,y - A5), 7 = 1....M (7) 

5=1 

The numbers are even or odd half integers depending on iV — M — 1 being 
even or odd, and together with the integers nj they specify a solution of 
eqs(3,7). 

Each allowed choice (see below) of {tcj, I^} uniquely determines an eigen- 
state of the Hamiltonian. We shall refer to the {rij^I^} configuration as the 
quantum numbers of the state they determine. 

These quantum numbers replace, for example, the {n'j,nj} quantum 
numbers of the free electron gas, where in the conventional Fock-basis each 
level rij could be populated by a spin up and a spin down particle. When 
the impurity is removed the Bethe ansatz equations describe a free electron 
gas. Indeed, in the absence of the term 0(2A) in eq(7) we have, 

AT^ 2 M 

E = EY^i + ^E[0(2A7-2)-7r] 
i=i 7=1 

= Et^. + 7E[-2^A + E0(A7-M] 

i=i ^ ^7=1 5=1 

i=i 7=1 

namely, a noninteracting gas given in a basis which is already charge-spin 
decoupled, and therefore adapted to turning on the spin exchange interaction 
which couples to the spin sector only. 

What restrictions are there on the choice of configurations {nj^I^} ? 
Obviously, the spectrum is unbound from below, as the integers nj can take 
arbitrarily large and negative values. To define the model we introduce a 
"bottom to the cutoff K , taken to be very large compared with any 

physical parameter in the theory. Then, since we are interested in the low 
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energy properties of the model, we may study it in the Hmit where the cutoff 
is taken to infinity as long as the physical quantities of the model have a well 
defined limit. We impose the cutoff as follows: 

I Y^il < (8) 

The cutoff K is imposed on the eigenstates of the fully interacting Hamil- 
tonian and thus differs from the conventional cutoffs that are imposed on the 
eigenfunctions of the free Hamiltonian. The choice of cutoff is irrelevant in 
models, such as the Kondo model, that are renormalizable. It affects the way 
physical scales (such as the Kondo temperature) depend on bare parameters, 
but not the way physical functions depend on physical scales [1]. 

From the action of the spin flip operators B{A^) on \<jO > it is obvious that 
the state B{Ai)...B{Am)\<^ > vanishes if two of the A's coincide. Further, 
as |0(a;)| < TT, it is clear that the must satisfy the restriction 

r{N,M) = -{N -M - l)/2 <I^<{N-M- l)/2 = I+{N,M). (9) 

We shall call a configuration {/-y}, for which a solution exists with all 
A-y distinct, allowed. Counting all allowed configurations (using the string 
hypothesis, see below) one finds that there are, indeed, 2''^ configurations, as 
required by the dimensionality of the spin space. 

We turn now to the determination of the eigenstates beginninig with the 
ground state configuration 

The ground state. 

The state with the lowest energy is a spin singlet, Mq = iV/2, induced by 
a configuration of consecutive {/°}, 

I'+i = + 1 (10) 
with the /-y filling all the slots from /"*" to /~ 

I~<I-y<I^, = ±{N/2-l)/2. (11) 

There are N/2 slots from /"'" to /~ and all are occupied by the Mq = N/2 
spin quantum numbers I^. 

The charge quantum numbers {n^} are taken the minimum allowed by 
the cutoff. They are all distinct and run from —KL/2'k upwards. Setting 

= we have K = 2t^D, D = N'/L. 



36 



We shall be interested in solving the equations in the thermodynamic 
limit: N^L ^ oo, with D = / L held fixed (later we shall also take the 
scaling limit: K = ttD — )• oo to achieve universality). Therefore, rather than 
finding the actual solutions {A-y} we consider their density cr(A), describing 
the number of solutions in the interval (A, A + dA). In other words, 

^(A,) = l/(A,+i-A,). (12) 

When all A-solutions are real (which is the case of the ground state) eqs. 
(5) and (7) can be rewritten in terms of the A-density 

27r r 

E = J2 —rij + D dAa{A)[Q{2A - 2) - tt] (13) 
j=i ^ 

and 

iV^0(2A^ - 2) + 0(2A^) = J JAV(A')0(A^ - A') - 27r/^. (14) 

An equation for the A-density in the ground state, (To(A), is obtained by 
subtracting eq(14) written for A-y from that written for A-y+i and expanding 
in the difference A A = A-y+i — A-y, assumed to be of order 1/N. One then 
finds 

^„(A) = /(A) - j K{A- A')ao{A')dA', (15) 



where 



/(A) = 2^ 

TT 



+ 



c2 + 4(A - 1)2 c2 + 4A2 



= -T— ri^^2(A). (16) 

Here we used 0'(A) = —2c/{c^ + A^), and the fact that, for the ground state, 
/-y+i — /-y = 1 for all 7. We also define, for future use, 

1 77- 1 

@4x) = 0(^^) (17) 
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The solution of eq(15) by means of Fourier transform is immediate, 



cosh^(A 



+ 



cosh -A 



(18) 



The transformation properties of the ground state are found by calculat- 
ing Mo = J (To(A)dA = |iV in accordance with the consideration in eqs(9,ll) 
for finite N and M. The state has a Young tableau of two equal-length rows 
and is a SU(2) singlet. 

The ground-state energy is given by 



Eo 



J2 J^nj + D J dAao{A) [0(2A - 2 

i 

T{l + ic)T{^ 



TT 

2L 



{N' 



iDln 



Til-ic)Ta + i 



tc] 



(19) 



To show that, indeed, this is the lowest energy state we shall study vari- 
ations from the ground-state configuration These correspond to 
excited states. 



Elementary Excitations 

Charge excitations (particle-hole) — obtained by exciting the charge de- 
grees of freedom. Thus excite a given where —K < {2Tr / L)n^ < 0, to 
n'j = -\- Atc > 0. The change in energy involved is 

27r 

AE = --An > 0. (20) 
Ij 

Obviously M, which depends only on the {/-y} quantum numbers, does 
not change and neither does the spin. 

We see that the charge spectrum of the theory is that of a decoupled free 
gas, a result of the interaction acting only on the spin degrees of freedom. 

Spin excitations — obtained by varying the {/°} sequence from its ground- 
state configuration, leaving the charge quantum numbers, tcj, unchanged. 
One way to modify the sequence is to put "holes" into it, namely to have 
unfilled slots, and correpondingly omit A's. Another way is to add complex 
conjugate pairs of A's. 

The triplet: The simplest excitation (keeping the number of electrons 
fixed) is obtained by considering the state with M = Mo — 1. This is a spin 
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triplet since S = N/2 — M = 1. Equivalently, AM = — 1 means that one 
box is moved from the lower to the upper row in the Young tableau. 

This choice induces two holes since now = ibiV/4 yielding iV/2 + 1 
slots for the M = N/2 — 1 /-y's leaving two slots unfilled. To find the effect 
of a hole suppose we choose a sequence {/-y} such that I-y^+i = I-y^ + 2, and 
/-y+i = /-y + 1,7 7^ 7o, omitting the integer J'* = I-y^+i- The spin momentum 
corresponding to it, A'*, constitutes a "hole". This means that we have to 
solve eq(14) in the presence of a (bare) hole density (t'*(A) = ^(A — A'*). 
To be more precise consider an {/-y} sequence with holes in it, denoting 
the omitted integers by {Ij ,j = l,...,m}. The BAE with the prescribed 
quantum numbers {/-y} determine the corresponding set {A-y}. Now define 
the function 

1 M 

uiA) = - — [iV^0(2A - 2) + 0(2A) - E ©(^ " A7)], (21) 

constructed with the determined values A-y. Those values of the variable A 
that satisfy 

u{A^) = I, , (22) 

where /-y is an integer belonging to the sequence, are the solutions we began 
with, while those values of A satisfying 

v{A)) = , (23) 

where the integers omitted from the /-y sequence, are the holes. Introduc- 
ing the distribution functions «t(A) and a (A) of the A-solutions and A-holes, 
respectively, we have 

^ = ^(A) + ^'^(A), (24) 

since the number of holes and A's in the interval dA is given, on the one 
hand, by [c(A) + (t'*(A)]c?A and, on the other hand, by the number of values 
of Ij and /-y which u(A) takes as it ranges over the interval dA. 

The equation for the density cr{A) in the presence of holes is obtained 
from eq(21) by taking the derivative with respect to A, 

(7(A) + (t''(A) = /(A) - J K(A - A')(T(A')dA' , (25) 
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where 



a\A) = j:SiA-A';). (26) 
i=i 



The solution (in Fourier space) to the equation is given by 

Hp) = ^o{p) + aHp) (27) 

with 



Aaip) = -y e-^^'^ ^r^' (28) 
^ 2coshfp ^ ' 



^AS exp-|p| 
i=i 2 cosh f J 

being the change induced in the density. 

Thus the "bare" hole density e~^^ ^ (in Fourier space) is "dressed" to 
g-»A''p ^exp^^|p|^ by the back flow of the sea of spin momenta A: since all A 

momenta are coupled through eq(14), removing one of them affects all and 
leads to the redistribution given by A(t(A). 

Given the density we can calculate the properties of the state; M, the 
number of "down spins", or the length of the lower row in the Young tableau, 
is 

( 11 
M = \ a(A)dA = a(p = 0) = -N - -m (29) 

so that each hole contributes (AM)h = —\ and corresponds to a spin- 1/2 
object, a spinon in modern parlance, since it obviously carries no charge. 

The triplet excitation, being characterized by AM = —1, is made of two 
holes, again in accord with the finite N considerations. In the language 
of spin representations the state consists of a symmetrized product of two 
spin- 1/2 objects yielding spin-one. 

The triplet excitation energy /S.E*' for holes at Aj* and A^ is given by 

AE' = D J A(t(A)[@(2A - 2) - Tr]dA 

= 2D tan-^ e(-/c)(Af-i) ^ ^D tan"^ e(-/c)(Aj-i)_ ^gg^ 

It is a sum of two terms, each term being the energy carried by the spin-half 
spinon. 

We claim, then, that the triplet is composed of two spin- 1/2 uncharged 
objects whose spins are coupled symmetrically to form a spin-1 state. To 
confirm this interpretation we need to show that another state exists where 
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the spins are coupled antisymmetrically to form a spin singlet. We are lead 
to consider configurations of inducing complex A-y solutions. Since the 
energy is real these complex solutions occur in conjugate pairs. We shall 
show that a singlet excited state is composed of a "sea" of real A-solutions 
with two holes, at and Aj, and a two-string^ namely, a pair of complex 
A's located at A=^ = A ± zc/2, where K = \{K^ + A^). 

The singlet: The equations governing the state are obtained from eq(7) 
written first for real A's (1-strings) and then for the A's in the 2-string. The 
first equation determines the density of real A-solutions «t(A) 

a{K) + a^{K) = /(A) - J K(A - A')(T(A')dA' - (t''(A), (31) 

where 

(t''(A) = S(A- A^) + S(A- A'^), 
a'\A) = Ks(A-A) + K^(A-A). 

The second equation 

M M 

N^Q{A - 1) + 0(A) = -27r/(2) + ^ 0,(A - A5) + ^ 03(A - As) (32) 

5=1 5=1 

fixes the position of the 2-string A. 

As before, (t'*(A) arises by placing holes in the ground-state (consecutive) 
sequence. The string term (T**(A) is the contribution 0(A-y — A~'~) + 0(A-y — A~) 
of the two-string at A^ = A it ic/2 to the sum in eq (7). This contribution 
can be rewritten as 0i(A — A) + 03(A — A) with A real, recall @nix) = @{^x). 
When we convert the set of algebraic equations to an integral equation, by 
the method described before, we obtain eqs(31). 

Its solution in Fourier space is 

Hp) = Mp) + ^^'\p) + ^Ap), (33) 

with 

^Ap) 
Aa^\p) 



exp^lPl 
2 cosh 



_g-(c/2)|p|g-iAp_ 



+ e" 
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When o-(A) is fed into eq(32) we find A = |(A5* + A2), as stated. 
This state is indeed a singlet, as can be deduced from calculating 

M = 2 + J dAa{A) = 2 + a{0) 

= 2 + ao(o) + Aa\o) + Aa'\o) = ^n, 

leading to a Young tableau with two rows of equal length. 

The singlet excitation energy, A£^*, is also easy to calculate, and is found 
to be degenerate (in the thermodynamic limit) with the triplet excitation 
energy, 

AE' = D J [A(t''(A) + A(t'*(A)) (@(2A-2) -Tr)dA 

+D (0(2A+ - 2) + 0(2A- - 2) - 27r) 

= 2D(tan-^e('^/=)(^f-^)+tan-^e('^/=)(^'-^)), (34) 

since / JAA(7^*(A)[0(2A - 2) - tt] + 0(2A+ - 2) + 0(2A- - 2) = 0. We find 
that the string induces a change in the sea of real A momenta which exactly 
cancels its direct contribution. 

We thus confirmed the picture of spin-1/2 objects, each corresponding to 
a "hole" in the sequence, being the underlying fundamental excitations in 
the spin sector [1]. But though the energy of a two spinon state does not 
depend on the total spin, it does not follow that the interaction between the 
spinous is spin-independent. Indeed, we shall calculate the phase shift when 
these excitations cross each other and find that they differ in the singlet and 
the triplet state. 

The doublet: To have a single spinon excitation it is necessary to add 
an electron to the system whose total spin will then change to be one-half. 
This change induces a hole in the sequence since M is unchanged and N 
is increased by one so that there is another slot available which is unfilled. 
Denoting by A'* the corresponding spin momentum we find that the energy 
associated with the adding an electron to form a doublet is composed of two 
terms 

AE'' = 2D tan-i e(-/c)(A''-i) ^ (35^ 

the first term being the spinon energy and the second the holon energy (n 
is the level into which the electron was inserted). In other words, the added 
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electron decomposes into two independent excitations, one carrying the spin 
content and one carrying the charge content of the bare electron. 

Another way of creating a one-hole state is by removing a particle from 
the system. Thus consider N ^ N — 1 and M — )• M — 1. The counting 
argument indicates that a hole opens in the sequence leading to a spinon 
excitation, but since an n level is removed [n is filled in the ground state, 
< ) we have an antiholon, a massless anticharged particle, rather than a 
holon 

AE'' = 2D tan-i e(-/c)(A''-i) + (35) 

Thus the states obtained by adding an electron have the identical spin- 
content as the states obtained by removing an electron. The decoupled charge 
sector is, of course, different. 

We proceed to show that this basic excitation, the spin-| spinon is a 
massless relativistic right-mover, and calculate its scattering phase shift off 
the impurity as well as off another spinon. From eq(35), we learn that the 
Kondo system has no spin-gap, since there are arbitrarily low-lying spin- 
flip excitations for A'* arbitrarily large and negative. The spinon excitation 
energy takes the form e = 2Toe^'^^°'^^ with To, given by To = De~'^^°. This is 
a dimensionful scale, dynamically generated by the model, and in terms of 
which all other dimensionful scales are measured. 

Having identified a physical scale. To, we may take the scaling limit, 
D — )• 00 with the scale held fixed. The coupling acquires a dependence on 
the running cut-off c = c[D) = tt/ In ^ — )• as D — )• 00. In this limit, having 
removed the cut off dependence, all surviving quantities are universal. In 
particular, the excitation energy is 

e = 2roe^ (37) 

with 9 = (7r/c)A'* held finite in the limit. The form of the energy is typical 
of a relativistically right moving particle with 6 being its rapidity. As such 
it is also its momentum, p = e. 

Within the scaling regime we still distinguish between (universal) ultra- 
violet behaviour as T, /i, e, ... ^ To and (universal) infrared behaviour as 
r, /i, e, ... <^ To- Of course, all energy scales are small compared with the cut 
off D which has been taken to infinity. We shall find soon that the limit- 
ing procedure assures a universal limit to all excitations and thermodynamic 
quantities. 
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We now wish to compute the various spinon scattering phase shifts. To 
deduce them we return to the cut-ofF theory on a large but finite ring X, 
and note that A'* is not a continuous, independent variable. It is determined 
by the {/-y} configuration, and takes values which are determined by the 
choice of the omitted J'*. This allows the calculation of the S-matrix, or 
the scattering phase shift S{6)^ of the spinon off the inpurity, as well as the 
scattering of one spinon off another, — 62). The method we use [8] is 
based on the well known observation that when a particle is quantised on a 
finite ring of length L, the shift of its momentum from the free value is 
interpreted as a phase shift, 

« 27r 8(9) TT . 

p = 2loe = —n -\ — , = —A (38) 

L L c 

We shall refer to this method of calculation as the momentum shift method 
and will use it also in the next lecture to calculate the S-matrix of the Hub- 
bard model. Another method of calculation is due to Korepin [9]. More 
recently a bootstrap Kondo S-matrix approach was developed by Fendley 
[10]. 

We determine the phase-shift as follows: consider again eq(21) 

u{A) = - — [N'@{2A - 2) + 0(2A) - / dA'@{A - AXA')] (39) 

and evaluate u(A) for a density (t(A) in the presence of holes {A^}, eq(27). 
We shall add 2-strings later. One finds 

27rz/(A)= - 2iV^(tan-^e^(^-^) --)-2(tan-^e?^--) 

2 2 

+ (40) 



Consider first the case with one "hole" in the {/-y} sequence, a one spinon 
excitation (induced, for example, by adding an electron to the system). The 
corresponding A'* is given by the value satisfying 27ri/(A'*) = j'*, where j'* is 
the "hole", so that eq(40) takes the form (upon dividing hj L ): 

2^(Un-' e?(^''-^) - -) = -/'^ + 2-(- - tan"^ e?^") (41) 

2' L L^2 I y I 
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In the scaling limit the left hand side contains the spinon momentum, 
and we find that its shift from the free value ^l'^-, namely the the scattering 
phase shift, is 



8{e) = 2(- - tan-^ e^) (42) 
and the spinon-inpurity S-matrix becomes 

^---1^— M^^). (43, 

The phase shift we found describes the one- dimensional scattering pro- 
cess. If we want to interpret it as an s-wave scattering in three dimensions, 
we have 

S = e''*" (44) 
with ^0 being the s-wave phase shift so that 

^o(p) = ^%) = |-tan-^(^). (45) 

As a function of the momentum the phase shift varies from = f at 
low momenta to = at large momenta. We shall interpret the the large 
momentum physics as being described by a weak coupling hamiltonian, iJo, 
while the physics in the infra-red is given by another hamiltonian H* which is 
strongly coupled producing a maximal phase-shift. This is our first example 
of a crossover in the properties of the model. 

The phase shift describes the scattering of the spin content of the incom- 
ing electron. The charge does not scatter. Thus, after passing the impurity 
the emerging state is not an electron anymore, but a superposition of un- 
scattered holons and a time delayed spinon. 

The spinon-spinon S-matrix can be read off eq(40) by considering the 
configuration with two "holes" at Aj* and Aj corresponding to the omitted 
integers and In this state the two spinon scatter in a triplet state. We 
have 

'"/'^-. +6(6)+ 1 m-^^)Til + ^^) 
L tL r(i + iAi_^)r(|-i^) 

with a similar equation determining p2. Hence, the spinon momentum, pi = 
Toe^i, is shifted by ^(^i) in the presence of the impurity, and by (f)*(9i — 62), 
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the triplet spinon-spinon phase shift, in the presence of another spinon. The 
triplet scattering is 



r(-{ _ jii^\r(l A_ 

s\e, - e,) = e^^'(^-^^) = i , i ^ (47) 



r(i + i^)r(| 



t- 



27r / \ 2 27r 



In a simialr way, to determine the singlet scattering phase-shift we have to 
evaluate eq(21) in the presence of two holes and a 2-string, using the density 
in eq(33). One finds 

S'ie,-e,) = ]~[^ s\e,-e,), (48) 

TT 

and hence the total spinon S-matrix in the Kondo model is 



sf. - - ^r(i+,:^)r(i-,^)j ^ + .0 ' 

The S-matrix we obtained again satisfies the YBE as does its bare counter- 
part. This guarantees that the spinon quasi-particles, although interacting, 
do not decay. In fact they are protected by the conservation laws. The same 
S-matrix was previously calculated for the spinous in the g-ology model [8], 
and will appear again in Lecture 4 to describe the scattering of the Hubbard 
spinous. The same S-matrix will describe the scattering of spinous in all 
models solved with the R-matrix of Lecture 2. It can be also deduced from 
general bootstrap considerations [11]. 

The spinon-holon S-matrix, S"''^ = 1, linearizing the spectrum has com- 
pletely decoupled the charge from the spinon sector. 

The full spin spectrum of the model can be constructed this way; excita- 
tions are built of holes and complex pairs. If we have only real A solutions 
with m holes in the A-sea, then the state will have spin S = |m. (With N 
fixed m must be even.) These are maximum-spin excitations. The energy 
associated with holes at Aj*, • • • , A'^ is 

m 

(AE)"" holes = 2^Dtan-i(e("/=)(^'-')). (50) 
i=i 
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When complex pairs are added they lower the total spin by coupling spins 
antisymmetrically. The various complex structures allowed will be discussed 
below, but they have the feature that their contribution to the energy can- 
cels. The latter therefore is determined only by the holes in the sea of real 
A momenta. The contribution of the complexes shows up, however, not only 
in the counting of states and their total spin but also in the S matrix deter- 
mining the interaction of the various excitations. We conclude that the spin 
excitations form a Fermi liquid. 

The nature of the complex solutions of the equations is captured by the 
string hypothesis: the solutions of eq(7), in the limit N — )• oo, always occur in 
the form of ra-strings, where an ra-string is a complex of n A-solutions given 

by 

Af = A('^)+i^(ri + l-2i), i = l,2,...,ri. (51) 

This hypothesis has a long history [15], [12], but is not always valid; the hy- 
pothesis holds in the presence of a macroscopic number of holes, and is not 
necessarily true in case we consider excitations containing only a small num- 
ber of them [13], [14]. Two-string solutions always exist (when two holes 
or more are present), but the conjugate pairs organize themselves into n- 
strings, > 2, only if driven macroscopically. Thus more care has to be 
taken when analyzing scattering events of elementary excitations than in 
thermodynamic applications, where, in studying the response of the system 
to external probes, we excite a sufficiently large number of holes and the 
string hypothesis is valid. 

Let us develop the form that eq(7) takes in this case. Consider the case 
with Mn TC-strings A^"j = A(") +z(c/2)(tc + 1 - 2j), 7 = 1,.. . ,M„, with A(") 
the real part of the 7th ra-string. Equation (7) then becomes: 

iV^04AW - 1) + 04 AW) = ^ 0.,™(AW - aJ™)) - 27r/W , (52) 

where the summation is over all strings different from the particular A^"^ 
string. The function Qmn{x) is defined as 

Q\n-m\(x) + 20|„_„| + 2(a3) H h 2Qn+m-2(x) + Qn+m(x),n ^ 171 

202(a;) + --- + 202n-2(a;) + 02n(a;), n = m 

(53) 
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and as a reminder, Q{x) = —2 tan ^(x/c) = 0ii(a;). Equation (52) is ob- 
tained from eq(7) by summing over all members of a string. For example, 

1 " r - ?7A - A^"h 
i=i S-=i c + t(A-A)>) 

n 

^0(2Af-2) = 04A('^)-1). 
i=i 

The integers /^"^ determine the allowed string solutions A^^j, and are the 
spin quantum numbers of the system. 

The rest is just as before. For a chosen set of configurations /^"^ one 

determines the corresponding A^^j complex spin solutions, and then forms 
the function 



ZTT 



iV^04A - 1) + 04A) - J2 0n™(A - A 



(54) 



with the A^™^ determined earlier, so that solutions of UniA^"^) = /^"^ are the 
allowed strings, while solutions, Aj^"\ of 

u4Af-^) = (55) 

where /^^"^ are the integers omitted in the sequence, correspond to ra-string 
holes. In the limit of iV — )• oo we may introduce ra-string density (Tn{A) and 
TC-string hole density (T^(A), which obviously satisfy 



dA 



<(A) + ^4A). (56) 



Combining the last expression with the derivative of eq(54) we find that 
the string densities obey the following set of equations 

oo 

/4A) = ^^(A) + J2 AnraTraiA), 
m=l 

where 

/4A) = N''K4A-1) + K4A), 

Anm = [\n-m\] + 2[\n-m\ + 2] + --- + 2[n + m-2] + [n + m]l57) 
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and [n] is the operator given by 



M/(A) = J K4A - A')f{A')dA 



(58) 



In terms of the string variables the energy can be expressed as 



E = DJ2 f dAa^A) [Q^A - 1) - ^r] + ^ 



(59) 



i 



Here we have performed the sum over the individual members of a string 
and are left with integration over the string locations only. 

With this complete classification of states, we proceed to compute the 
partition function of the system at nonzero temperature T and external mag- 
netic field h. We shall deduce a set of coupled integral equations determining 
the free energy using a method originated by Yand and Yang [16], and 
generalized by Gaudin [17] and Takahashi [12]. We shall analyse the equa- 
tions and demonstrate scaling behaviour and crossover properties in the full 
h — T plane. 

The Thermodynamics of the Kendo Model. 

The formal expression for the partition function is 



where H is the zero-field Hamiltonian, fj, is the magnetic moment and we set 
g = 2 for convenience. Tvgs^ is the trace over all basis states with values S 
and Sz of total spin and z component of the spin. Since H is invariant under 
simultaneous rotations of all spins, we may split off the sum over to obtain 



Z 





E E 



(60) 



5=0 5^ = -5 



z 



iv/2sinh {2S + l)f 

E f — 1 — 

s=o sinh ^ 



Hi 



T 



N/2 

J2 Tvss exp 




(61) 



5=0 
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In the last approximation we have dropped terms proportional to exp( — S'/i/T), 
as well as an overall factor [2 sinh(/i/r)]~-^, since these terms contribute neg- 
ligibly to the calculation of thermodynamic quantities in the limit X — )• oo 
(note that S~L). 

Now let us exploit the specific form of the energy for our basis states. 
Recall that each such state is labeled by a set of quantum numbers {rij^I^} 
with rij > —N^. The corresponding energy is 



with the charge energy 



i=i 



(62) 



(63) 



and the spin energy 

£;W({A}) = dy: 



M 



7=1 



0(2A, 



TT 



The partition function factorizes acordingly, 

where the charge partition function 

= J2 exp 

{nj},nj>-N- 



T L \ 



(64) 



(65) 



(66) 



describes the thermodynamics of N'^ noninteracting spinless fermions with 
linear kinetic energy. In the limit D — )• oo it leads to the free energy 



= -— / dk\n{l + 

27r J-oo 



_k_. 

e T 



TT 



— + {infinite constant} (67) 



and is half the the free energy of a noninteracting electron gas at zero mag- 
netic field. Note that the effects of a magnetic field are all included in Z^*). 
The spin partition function is 



exp 



V T J 



I] I] exp 

M {Ai,...,Am} 



-[£;W({A}) + 2MM] 



(68) 
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It can be rewritten in terms of the ra-string and ra-string hole density «Tn(A) 
and cr^(A) as 



Z(') = exp 
with 

where 



N/ih 



J n^o-„^o-^exp<S({o-„,o-^}) exp 



E^'\{A}) + 2fihM = J2j dAa4A)g4A) 



[E^'\{A})+2nhM] 
(69) 

(70) 



gn(A) = D[Q4A - 1) - tt] + 2fihn (71) 

and where iS({«t„, «t^}) is the entropy associated with the densities {cn,cr^}- 
In other words, exp iS({«t„, «t^}) is the functional, counting the number of 
configurations {/-y} leading to solutions {A-y} that are consistent with a given 
set of densities {cr^, To determine S divide the A axis into intervals dA, 
chosen sufficiently small so that the densities are approximately constant 
over each interval, yet sufficiently large that (cr^ + a^)dA ^ 1. The number 
of slots for A-strings in the interval dA is dun = («Tn + fn)'^A, a^dA of which 
are occupied, while cr'^dA are empty; thus the number of ways of distributing 
the TC-strings among the slots is 



[(^(A) + ^^(A))JA]! 
[a^{A)dA]\[ah{A)dA]V 

Using Stirling's formula, we can simplify this to give 

[{a^ + a^JdA]\ 



(72) 



dSr, = In 



[a^dA]\[alldA]\ 
so that entropy, iS, becomes 

S({<rn, <rt}) = E / d^i^^n + <rt) Han + a 



[{ar^ + (7^) ln((7„ + (7^) - (7„ lu (7„ - (7^ lu (7^] JA 



C7. 



In 



C7„ 



In 



(73) 



In thermodynamic limit, iV — )• oo, we may evaluate Z^*) by the method 
of stationary phase approximation. Varying the functional 



^^'H^n, ^T^} = E^''^ + 2nhM -TS-Nh 



dA 



a^g-a - Tan In 



1 + 



C7' 



Ta'^ln 



1 + 



a' 



(74) 
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subject to the constraint, ^cr^ = — Anm^f^m-, one finds, 

fl'n(A) 

-± 

m=l 

where 



+ ^ A_ln[l + 7/-i(A)], 



(Tn(A) 

This set of equations may be rewritten as 

ln7/„ = G'[ ln(l + 7/„+i) + ln(l + 7/„_i)], 
on 

In 7/1 = -— tan-'e("/=)(^-') + G'ln(l + 7/2) 



(75) 
(76) 
(77) 



with the driving term, 2D tan ^ e^'^/")^''^ ^\ familiar as the one-hole excitation 
energy AE"^. The operator G is defined by 

[1] 



GfiA) 



[0] + [2] 



cosh^(A - A 



7Y/(A')- (78) 



To close the set of eqs(77) one has to supply a boundary condition for 
— )• oo. It turns out to be 



lim {[n + 1] ln(l + 7/„) - H ln(l + 7/„+i)) 



2iJ,h 



(79) 



Once a set of solutions r/n satisfying the equations has been found, the 
spin free energy may be obtained from eq(74), and considerably simplified 
(we shall skip some of the steps). 



nm 0"m ln(l + 7/„^] 



a^Qn - ro-„in(i + ?/„)- r 

Er / JA/4A)ln(l + 7/-^)-iVM 

J2T f dA[ N'S(1 - A) + S(A)][n]ln(l + 7]-^) - Nfih 



ln(l + 7/i] 



-r[iV^^(l- A) + ^(A)]G' 
I dAao{A){g^{A) - T\n[l + r]^{A)]} - Nfih 
Eo-T I JA^7„(A)ln[l + 7/i(A)], 



N/ih 



(80) 
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where 



a^iA) = G[N'Sil - A) + ^(A)] = ^ 

Zc 



+ 



cosh^(A-l) cosh^Aj 



(81) 



is the ground-state A density, and Eg is the is the ground state energy. 
Adding the charge free energy FW, eq(67), we have 



F = + F(^) = Eo- -T j JA^o(A)ln[l + 7/1 (A)], (82) 

Eo now containing also the temperature-independent contribution of the 
charge fluctuation. 

We shall now demonstrate the scaling properties of the thermodynamic 
equations, describing the behaviour of the system in the regime where T^h « 
D. In this case the function 7/1 has a very sharp decrease, proportional to 
exp[— (2Z)/r) tan"-^ z], where z = exp[(7r/c)(A — 1)], and will contribute of 
order exp(— 2Z)/r) to the partition function except for small z. Hence in the 
scaling limit we may replace tan~^ z by z in these integrals and compute rjn 
from a modified version of the thermodynamic equations 

ln7/„ = -2^„,ie'^ + G'[ln(l + 7/„+i)+ln(l + 7/„_i)] r^ = l,2, (83) 

where we now regard rjn as a function of the new variable C, 

C = ^A-lnf, (84) 

with 

and by convention we set 7/0 = 0. 
The free energy becomes, 

ttXT^ T /• , f 1 1 , , ^h^^ 

^ = ^'--^-^ / '^^ 1 .1 + 1 r 

12 27r J y cosh ( - In ^ - ^ cosh ( - In ^ J J- 

(86) 

which we rewrite now in terms of impurity and electron contributions rather 
than in terms of spin and charge sectors, 

F = Eo + F' + F\ (87) 
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The impurity free energy 



depends on the coupHng constant and the cutoff only through the combina- 
tion To = Z)e~'^/°, the only scale in the problem. The free energy is universal 
to all materials or constructions with the same Tq. 
The electron contribution 



12 




(89) 



is the free energy of a non interacting gas of electrons, here expressed in the 
charge-spin decoupled basis. This is obviously so since it is the part of the 
free energy that survives when the impurity is removed. 

Our task, then, is to solve the equations eq(83) and from the obtained 
Vi{C) to compute the impurity free energy F*, eq(88). Some properties of 
the solutions are easy to establish, even though no analytic solution is yet 
available: 

1. 7/n(C) is monotonically decreasing in ^ ( fixed n). 

2. 7/n(C) is monotonically increasing in n (fixed (). 

3. i]n{() has finite asymptotic limits: 



I sinh? (n+l) 'if!- ^ ^ 

Vn = — ■ .2 M - 1, as C ^ +00 

y (90) 

_ sinh^ n ^rr ^ u. ^ ' 

= - 1' as C ^ -oo 

These properties will be of use as we turn now to study the impurity 
physics in the high tempertature and low temperature limits for a given 
magnetic field and then at zero temperature as a function of h. 

Let us check that the electron free energy F'^ coincides in the limit D — )• oo 
with the conventional expression derived in the Fock basis [2]. Begin from 
the expression for F^""^ = — J2n T J dAfn{A) ln(l + t/""^), and notice that now, 
as we are considering a free electron gas does not contain the impurity 
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part, /n(A) = N^Kn{A — 1), which is proportional to the derivative of g„ 



/n(A) 



' 2ir 



(/^(A) which in turn via eq(75) allows us to write 



E T'f^ J dkj^ [ln(l + 7/.) - E ln(l + 7/-^)] ln(l + t/"^] 

y.T'- [ 



ln(l + 7/„)-ln(l + 7/-^) ln(l + 7/-^] 



-ET^^ /'^A^ln(l + ^-^) = -Er'^ J,.^ln(l + ,-) 

X 1 X /''^ 

-T^— dr]-\n(l + r]-^) = -T— dk\n(l + e'T). (91) 

27r Jo J] 2'ir J-oo 



We used the fact that rj^ = 0, 7/^ = and in the last step changed 

k 

variables rj = ct ^ to obtain an expression recognizable as the free energy of 
a gas of non interacting spinless electron, and identical to F^°\ both being 
half the conventional free energy. We turn now to investigate the impurity 
physics. 

Impurity behavior at high temperature. 

The impurity free energy at high temperature is determined by the func- 
tion 7/1 (C) for large negative values of (. In the limit it is found from r/^ to 
be 



T 
'2 



TT J 



— — ^— ^In [l + r^,{C^^)] 
cosh [C + In T 



-Tin 



2 cosh 



fj,h 



(92) 



This is the free energy of an isolated spin in the presence of a magnetic 
field /i, and we are therefore in the neighborhood of the weak coupling fixed 
point. How rapidly is this point approached? Assuming an expansion of r/n 
in inverse powers of ( and In C, inserting the expansion in the thermodynamic 
equations we can determine the asymptotic behaviour of r/i and hence the 
impurity free energy 



ln(2 cosh 



fj,h . 



1 fj,h fj,h 
— -^ tanh — - 

2 T T 



l lnln(r/rfc) 
In T/Tfc ' 2 WTjTk 



+ 



(93) 



leading to the susceptibility. 
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^ TV Vln T/Tfc 2 WT/Tk ) V In" T/T^ . 



(94) 



We have introduced a new scale the Kondo temperature Tj,, defined by the 
requirement that no term of the form l/ln^(r/rfc) appear in the expansion. 
In other words, the term a/ln^(r/ro) that actually does appear is absorbed 
by expressing the expansion in terms of a scale = cTq. This expression for 
the susceptibility can also be obtained by conventional perturbation theory, 
by studying the neighborhood of the weak coupling hamiltonian, iJo, the 
ultra-violet fixed point. 

Impurity behavior at low temperature. 

We shall see now that as T tends to zero the system flows to a new 
fixed point that is Fermi liquid in character. Indeed, note that the function 
/(^, x) which in eq(89) describes for small t the physics of a free electron gas 
on all scales, also appears in eq(88), where it captures for small t the low 
temperature physics of the impurity. Since the former is the prime example 
of a Fermi-liquid it follows that so is the latter. More precisely, the impurity 
physics at low temperatures is the same as the spin sector of the free electron 
gas, and as such it carries all the spin degrees of freedom but only half of the 
entropy. This is the origin of the value of the low-temperature Wilson ratio 
R to which we turn. 

Let us then carry out the limiting procedure ^ = ^ ~^ 0. For sufficiently 
small t we may expand the kernel 1/ cosh(C + In ^) = 2t exp ^(1 — t^ exp 2^ + 
exp4C + •••), to obtain (we may perform the expansion within the integral 
as 7/1 vanishes for large ( as e"^'^'") 

fit, x) = - j dCei ln(l + 7/i(C, x)) + 0{t'), (95) 

TT J 

SO that 



L 12 



-/jCe'^ln(l + 7/i(C,-))+0(— ). (96) 



As argued earlier, is the free energy, in the spin-charge decoupled 
Bethe basis, of a system of = DL noninteracting spin-| electrons in the 
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presence of a field h. Comparing it with / L calculated in the conventional 
manner in the the Fock basis (the identification is valid as D — )• oo), 



F' _ T 



/OO fOO 
JA;ln(l + 6-*=/^) + / JA;ln(l + 6-*=/^] 

-(irD-fih) J-(irD-fj,h) 

/,Ce<l„(...,4,)=g.M£. 

The same integral appears in the impurity low-temperature expression, 
therefore we have 



we obtain 



F^ = -^JdCe^ln{l + ,.iC-)) 



—T' + -{^lhy 



7rroLl2 2 



(98) 



We find that the impurity contribution to specific heat at low temperature 
is Fermi-liquid like 

C = (99) 

as is the magnetic susceptibility 

= 4r- (100) 
Trio 

The zero temperature susceptibility is finite, indicating that the impurity 
spin manifest in the high temperature regime by Curie's law, = is 
now completely screened. We shall interprete the effect as due to strong 
effective coupling between the impurity and the conduction electrons leading 
to the formation of a local singlet, and the infra-red physics is dominated by 
a strong coupling fixed point. 

Although the strong coupling impurity physics is Fermi-liquid like, it 
is different from the electrons' Fermi-liquid. This can be brought out by 
comparing 



2 



^ = (101) 
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with the corresponding electronic value 
R = -^riTTT takes the value. 



3^. Hence Wilson's ratio 



R 



(102) 



The main element that allowed the identification of the fixed point as 
a Fermi liquid was that the same function r/i occurs in the description of 
both the electron gas and the impurity. At the same time we have R = 2 
due to the decoupling of the charge degrees of freedom from the impurity. 
More generally, the natural basis in Hilbert space for the description of the 
infra-red is the Bethe-basis with charge and spin decoupled, and the impurity 
modifying only the spin sector, while in the ultra-violet the Fock basis is the 
natural one. 

Let us study the crossover as a function of the magnetic field at zero tem- 
perature. In this case the infinite set of themodynamic equations collapses 
into a single one. 



/•oo 

aB(A) + K(A- A')aB(A')dA' = /(A), 
Jb 



(103) 



where the lower bound B is related to the magnetic field h. The magnetiza- 
tion equation can be derived directly by looking for the lowest energy state 
for a given spin. In this state the spinous are excited at the lower end up to 
A = B. 

The equation can be solved by means of Wiener-Hopf technique and yields 
the mag netization, M = fiS = ii\{N - 2M), as a function of the magnetic 
field [4] [3]. The solution is of form. 



M = M' + M' 



(104) 



where Ai^ = fi, 
magnetization o: 

2 11/2 



nl/2 



LTqC^^^^ in the scaling limit. Identifying it as the 
■ the free electron gas, MpauH = fitiL/ir , 

IJ^^Q^B/c ^ ^ ^ 2)) (105) 

LTreJ TT 



allows us to relate the parameter B to the magnetic field 



^ttB/c 



ire 



(106) 
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The second terms is the impurity magnetization which upon the previous 
identification becomes. 



Ti < nh « D. 

(107) 

This expression is vahd over the entire range of the energies and we may read 
off the asymptotics: 
In the ultra-violet. 



+ 



In 2 



In^ 

J-i 



21n^ In'^ 21n'^ 

J-i J-i J-i 



+ 



(108) 



which is the magnetization of a free spin weakly interacting with the conduc- 
tion band. This is in accord with our previous conclusions. 
In the infra-red[l], 



ttTo 



-h 



(109) 



indicating a screened impurity. 

This magnetic crossover also figures in another interesting quantity, the 
scattering phase shift of an electron on the Fermi surface off the impurity 
as a function of the magnetic field, ^o(^) [18]. This quantity determines the 
transport relaxation time r 



1 

27 ~ " 

and hence the magnetoresistance 



pih) 



-c sm 



3c 



sin^ So{h). 



(110) 



(111) 



Here c is the impurity concentration, and Dq the density of states at the 
Fermi level. 

To calculate ^o(^) we add an electron to the system which is at its lowest 
energy state in the presence of a magnetic field, lO/j >, creating a doublet, a 
one-"hole" state \k^a;i}h >■ Denoting by Eo{N^ ^ S) the lowest energy for a 
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system with electrons and a total spin S', we have that the energy of the 
state a; iih > with k on the Fermi-level is 

Eik, ±1/2; nh)\k=kHh) = E{N' + 1, 5' ± 1/2) (112) 

where S is the spin corresponding to the field h. On the other hand we may 
add the electron far from the impurity, hence E{k, it; O/j) = k-\- E{^h) hence 
only those values of k are allowed (on a finite ring L) which satisfy 

k = E(k, ±1/2; ilh) - Eiilh) = E(N' ± 1, 5 ± 1/2) - E(N% S) (113) 

By studying the deviation of the allowed incident momentum k from free 
values we can determine, as before, the scattering phase shift, 

So,±{h) = So,±{kF{h)) = ^[1 ± M\h)/fi]. (114) 

where the sign indicates whether the electron spin projection is parallel or 
antiparallel to S. 

The phase shift reaches its unitarity limit, = '?r/2, at low magnetic fields 
and falls to zero logarithmically as the fields are increased. The corresponding 
behavior in the resistivity (albeit as a function of the temperature) was the 
experimental measurement that ignited the interest in the model. 

We have explored the crossover as a function of temperature, magnetic 
field, and momentum concluding consistently that all infrared physics can be 
described by a strong coupling fixed point that produces a Fermi liquid. Let 
us turn now to discuss it further. 

The strong coupling fixed point. 

In the infra-red regime the system can be viewed as a gas of quasi parti- 
cles scattering off a non-magnetic potetial center characterized by the phase 
shift S{p)^ and weakly interacting among themselves. These are the spinous 
interacting via scattering matrix elements S'* and S'* and undergoing a phase 
shift So{p) upon passing the screened impurity potential. The localized po- 
tential is due to a singlet formed by the impurity strongly interacting with 
the electrons in the infra-red and effectively capturing an electron to form 
a local singlet. In other words, from Friedel's sum rule we can argue that 
So{p = 0) = 7r/2 indicates an enhanced density of states around the impu- 
rity, tantamount to an electron (or its spin content) captured there, forming 
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a singlet. Virtual transitions are then responsible for inducing the scattering 
of the spinon off the screened-impurity. 

We are immediately led to the following expession (valid at low temper- 
atures) for the impurity free energy 

„ I dp ^, , ^ , , P , TT 



' r^S'ip)lnil + e-^)^-^'- (115) 



as sum over spinous (hence the lower limit of the integral is zero). Here 
^^'[e) = = ~ (2To^°+e^ impurity part of the one-particle density of 

states. (An equivalent way to obtain it is from T> = ^ = (To(A)/ 

More generally, Nozieres[19] has argued that the physics around this fixed 
point can be captured by a Landau expansion of the phase shift 

^o(e) = 8Q-ae± (f)''2M^ (116) 

with ^0 = f from strong coupling arguments, and a and cjf' phenomenological 
parameters. These parameters determine all low temperature proprties of the 
model. In particular, resistivity specific heat and magnetic susceptibility are 
given by 



p{T) = po(l-(«7rr)' 
9i - ^ 

where Dq denotes the density of states, here Dq = ^. Wilson's Ratio is given 
by i? = 1 + 2Do(jf' I a. Assuming the Kondo singularity is tied to the Fermi 
level and only antiparallel spin interactions Nozieres concluded 2Do(jf' j a = 1 
and hence i? = 2 in accord with Wilson's result. 

The explicit values of these parameters can be read off eq(45) and eq(114). 



1 , TT TT TT 



2To 4:M^ ^TqL 

yielding a complete characterization of the fixed point. 
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The Kondo Problem. 

We now return to consider the main question: the characterization of 
the crossover by universal numbers [5] [20]. Begin by introducing the high 
magnetic field scale T/j, parametrizing the weak coupling regime in the h — 
T plane, by requiring that no term of 0'^'n~'^ {h j Th)] should appear in the 
expansion. One finds from eq(108), 



with 



In In 7^ 1 „ 



21n^ 41n'^ Un ^ 

J-h J-h J-h 



TH=(^yT,. (117) 



The number Wh = 3^ = y 7 [3] , defined as the ratio of the ultraviolet scale 
Th to the infrared scale To, relates to both fixed points and its computation 
requires the construction of the model on all energy scales. 

Can one similarly understand the crossover as a function of the temper- 
ature? One wishes now to calculate = where Tk was defined at the 
weak coupling regime and characterizes the corrections to Curie's law, while 
To defined at strong coupling sets the screening scale. We have encountered 
the number previously in the form W = e"' va. the discussion after eq(94), 
but were unable to determine it from the infinite set of coupled thermody- 
namic equations. We proceed to calculate it now, appealing to the idea of 
universality. Expressing W as 

H^=^ = ^^, (118) 
To ThTo' ' 

we note that the universal number Tk/T^ is completely defined in the weak 
coupling regime, and therefore can be calculated exactly by means of pertur- 
bation theory. We find, 

T, 

2/37v^e-4 (119) 



Tk 



T, 



h 



where 



ln(3 = / dx{l-x^)x {—^ ^ =ln [21] 

Jo ysin irx x^ / 27r 

In 7 = C = Euler's constant. 
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Hence 

— = = 0.102676... [31. (120) 

47r 47r 

This number was first calculated by RG techniques to be 
W 

— = 0.1032 ± 0.0005 [5], (121) 
47r 

in good agreement with the analytic result. Note that the three numbers 
TkITh, Tfi/To^ Tk/Tq were computed within three different constructions of 
the model. Their accord is due to the renormalizability of the model and the 
universality resulting from it. 

We have seen that the Kondo model flows toward an infra-red behaviour 
that is Fermi-liquid in character. We wish to describe now some generaliza- 
tions where this behaviour is modified. 



The multichannel Kondo model. 

The model was introduced by Nozieres and Blandin [22] to include the 
orbital structure of the impurity, 

H = -ij ra,raix)d.i^a,rai^)dx + 7^,^(0) ^afe^fe,™(0) • S (122) 

here m = 1, ...,/ = 2/ + 1 is the orbital channel (or flavor) index and the 
spin operator S'* is in spin-S' representation of SU(2). In the hamiltonian the 
values of / and S are unrestricted. To describe a magnetic impurity / = 25', 
but other non-magnetic applications of the model exist with other values of 
spin and flavor [23]. 

The nature of the infra-red fixed point depends on those values [22] : for 
/ < 25' the coupling J flow to infinity leading to a screened impurity in the 
case / = 25', and to a partially screened impurity 5'' = 5' — | in the case 
/ < 25'. The strong coupling fixed point becomes unstable when / > 25' and 
the infra-red physics is then controled by a new, finite coupling fixed point. 
This new fixed point is expected to describe non Fermi-liquid behavior. 

The model is integrable but presents some challenge as to the handling 
of the cut off scheme. The reason is that the effect of flavor enters through 
the Fermi statistics; it allows / spin-up electrons to interact at a time with 
a spin down impurity. The cut off procedure, if improperly handled, may 
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smear out the local interaction and lead to wrong results. At the same time 
it must respect integrability. 

One procedure [25], valid only for / = 25', uses a version of the Anderson 
model as a cut-off. A more gereral approach is not to linearize the spectrum 
ab initio, but to maintain some curvature which is removed at the same time 
as the cut-off. This way one obtains a set of Bethe- Ansatz equations coupling 
spin and flavor degrees of freedom [24], and a remarkable effect takes place 
when the cut-off is taken to infinity: the flavor singlet sector separates from 
the rest of the Hilbert space leading to a new "fused" set of BAE describing 
the interaction of effective spin-//2 electrons with the spin-S' impurity. This 
dynamical fusion to form the spin-//2 electron complexes, captures the full 
spin content of the model and underlies the appearance of the non-Fermi- 
liquid behavior. 

We shall skip further details of the Bethe- Ansatz construction and write 
down directly the BAE for the model. 



il A,-A, + ic _ / A,-l-i/c/2 \^ l A,-icS \ 
l}-^As-A^-ic \A^-l + ifc/2) \A^ + icS)' ^ ' 

and the resulting thermodynamic equations in the scaling limit, 

\nr]n = -2^„j + G'[ln(l + 7/^+1 )+ln(l + 7/„_i)], (124) 

noting that the effect of flavor is to move the driving term to the /th row 
of the coupled set of equations. The asymptotic conditions on the equations 
are, 

[n + l]ln{l + rJr^) -[n]ln{l + rJr^+l) ^2 fih/T, n^oo. (125) 
Having solved the equations one computes the spin-S' impurity free energy 

from 7/n=25, 

(120) 

^'^ 27r J ^ cosh [C - In ^] ^ ' 

The same set of equations determines the free energy for all values of the spin 
S. We shall see below that the overscreened solutions n < f differ in 
character from those with n> f . The electronic properties are still calculated 
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from r/f. Hence we may expect Fermi-liquid behaviour only in the case / = 
2S. 

The functions r/n are analytic monotonically decreasing in ( for all n and 
tending to finite limits 7/^ as ^ — )• ±00. The limits are given by 



sinh (fih/T) 
and 



,mh'[in + l)^lh/T] 

, TC = 1,2,..., (127) 



( sin^[(n+l)7r/(/ + 2)] _ for « < f 

Vn - sinh-[(n+l-/y/T1 _ ^ fo^ „ > / ^^"'^ 
I, sinh-' (^J,h/T) ' — •' 

Consider now the high-temperature properties of the model. They are 
determined, for a spin-S impurity, by behavior of 7/25 in the limit ( — )• —00. 
As in the one-flavor case we have discussed in detail thus far this limit is 
approached with power corrections leading to 

, sinh(2S' + l)iih/T , , 

Lh,/./r 

This is the weak-coupling regime, governed by the fixed point at J = 0. 
The free energy is that of an isolated spin S up to the usual logarithmic 
corections. The nature of this fixed point is unaffected by the introduction 
of flavor degrees of freedom. 

On the other hand flavor affects significantly the low temperature prop- 
erties of the model. These are determined by the behavior of of 7/25 in the 
limit C — )• +00. The nature of the limit and the approach to it depend on 
the flavor degrees of freedom. 

The underscreened functions, n = 2S > /, attain their asymptotic 
limit 1]^ with power corrections and we have, 

Tl ^inh(25 + 1 - /)M/r 

This is the free energy of a spin S' = S — |/. In other words, the impurity 
spin is partially screened. The approach to the limiting value is logarithmic 
in the temperature indicating a fixed point at J = 00. 
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The function n = describes complete screening and, as T — > 
0, ~ (DiT^ + Sh^)/To. yli,i?i,Ci, and Di are numerical constants. 

A new situation arises when we consider the O'uerscreenec? functions, n < 
f. These functions approach their limit values exponentially, 7/n(C) ~^ + 
c„e~'^^,C — )• oo,n < /, with < r < 1. Hence we may apply the opera- 
tor G directly to the asymptotic form G e""^^ = e~'^^/2 cos(7rr/2) and the 
thermodynamic equations reduce to an algebraic recursion relation for the 
coefficients = c„/(l + ?/^), 

(&n+i + &n-i) = A&„, r^ = l,...,/-l (131) 



1 + 7/, 

with the boundary conditions bo = bf = 0. The eigenvalue A is related to r 
via A = 2 cos ttt /2. 

The solution to the resursion relation is 

sin[(r^ + 2)7^/(/ + 2)] 

sm[r^V(/ + 2)] ^'''^ 

with the smallest eigenvalue 

r = ^ (133) 

(In the two channel case more care needs to be taken since r = 1 and applying 
G to the asymptotic form generates logarithmic corrections.) 
The impurity free energy at low temperature takes the form , 

4. = -ItMi + ,i) - T,r/r„)' + . . . (134) 

from which we conclude that the infra-red physics is dominated by a new 
non-trivial fixed point generating power laws behavior for the specific heat 
and susceptibility 

C ~ i^)^ (135) 
1 T ^ 

X' ~ |(^''' (136) 
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A similar power law arises when one considers low-energy scattering of the 
spinon off the impurity. 

These power laws allow the identification of the infra-red fixed point 
hamiltonian as the SU{2) level-/ WZW conformal field theory where the 
combination is the coefficient of the hamiltonian in the Sugawara-form[27]. 
This coefficient determines essentially the conformal dimensions, and hence 
the the asymptotic form of the dynamic correlation functions [26] [28]. In 
particular, the primary field $i(^) has the dimension A = 2/(/ + 2) and can 
be identified as the spin operator from eq(136). 

Actually, the appearance of some WZW model as the fixed point hamil- 
tonian can be expected from symmetry arguments since it is the simplest 
model to incorporate scale invariance and SU{2) symmetry. Detailed calcu- 
lations are required to determine to which level (here level /) the full model 
will flow to in the infra-red. 

The zero temperature entropy can immediately be read off eq(134). 



entropy is due to the solitonic nature of the excitations [10], and is general in 
conformal field theory [29]. 

The fixed point can also be approached as a function of the magnetic field 
at zero temperature. The thermodynamic equations collapse into a single 
equation which describing the maximum spin excitations above the ground 
state, consisting of an /-string configuration in the presence of a magnetic 
field. For small magnetic fields, however, we can deduce the asymptotic 
behavior from the ground state solution. Denoting by ctq the impurity con- 
tribution to the ground-state density of / strings, we have for the impurity 
magnetization fields h 




sin[(25 + l)7r/(/ + 2)] 
sin[7r/(/ + 2)] 



(137) 



and is the logarithm of a fractional number! The appearance of fractional 




(138) 



The density ctq is given in Fourier space by 




sinh(S'JP)/[2cosh(Jp/2)sinh(/Jp/2)], f > 2S 
exp[(//2 - S)J\p\]/[2coshiJp/2)l f < 2S. 
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In the limit /i — )• 0, the magnetization is dominated by the properties of 
at p = 0. While for / < 25', o"q(p) is discontinuous at p = leading to 
-^(/i~o) ~ 2//(S' — |/) + O(lnh/To), for / > 25' the transform is analytic in p 
so that is controled by the pole at p = —2i/f. Hence 

M'(h) ~ n(h/ToY^^, h^O (139) 

leading to the critical exponent S = f/2 (and logarithmic corrections for 
/ = 2). 

We have dealt thus far with the theoretical aspects of the Kondo model. 
It is, however, a model that describes an experimentally realizable system, 
and the mathematical structure we discussed can be confronted with reality. 
It is found that the theory provides a remarkably good fit for a large body 
of experimental data [30] with no adjustable parameter except the scale Tk- 
This is more remarkable still in view of the simplicity of the model, which 
nevertheless captures the essential physics of a rather complex system. 
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Lecture 4: written with Andres Jerez 

The Hubbard Model 

In this lecture we analyse the spectrum of the Hubbard model. Many of 
the results have been derived in the past by Lieb and Wu [1], Ovchinikov [2], 
Coll [3], Shiba [4], Takahashi [5], Woynarovich [6], and Kluemper et al.[7]. 
We systematize the results, occasionally correct them, and derive new ones 
[8]. The latter category includes results concerning the spectrum and scat- 
tering of the excitation away from half filling, for the repulsive and attractive 
interaction. While these lecture notes were written up Korepin and Esseler 
also derived some of the results presented here [9]. Their work concentrates 
on the half filled case and emphasizes the charge SU{2) symmetry that ap- 
pears in this case [10]. This symmerty is explicitly broken away from half 
filling but we shall find some remnants of it. 

The analysis in this lecture is similar to the one carried out in Lecture 
3, but is much more involved due to the coupling of the Bethe-Ansatz equa- 
tions. We find that the fundamental excitations are, as in Lecture 3, spin-| 
uncharged spinons and spinless charged holons and antiholons. Indeed, any 
SU(2) invariant integrable model characterized by the R-matrix discussed 
in Lecture 2 will have excitations carrying the same spin quantum num- 
bers, since these are determined by counting. Also the spinon S-matrix, 
computed from a counting argument, has essentially the same form as we 
encountered earlier. The dynamics, on the other hand, captured by the 
function = determines the dispersion of the excitations and de- 

pends on the particular form of the interaction. In the repulsive Hubbard 
model we shall find that spin excitations are always gapless and have a spin 
Fermi momentum kp = ^n, where n = N/ L is the particle density. The 
charge sector is composed of gapless (holon-antiholon) as well as of gapful 
excitations (holon-holon). The gapless charge excitations are present only as 
long as the band is less than half filled, and their charge Fermi momentum is 
kp = irn. When n = 1 no gapless charge excitations survive and the system 
ungergoes a metal-insulator transition. 

In the attractive model the situation is somewhat reversed. Gapless 
charge excitations always exist: the holon-antiholon excitation is gapless but 
disappears from the spectrum at half filling, while the holon-holon excitation 
becomes gapless at this point. On the other hand the spin excitations are 
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always gapful. 

The nature of the spin gap can be easily understood by comparing eq(44) 
and eq(46) in Lecture 2. The couplings of the (antiferromagnetic) Kondo 
model and the attractive Hubbard are of the same sign. But as we saw in 
Lecture 3, the effective spin coupling then flows to the strong coupling fixed 
point generating a spin scale for the Kondo model, and a spin gap for the 
attractive Hubbard model. In the repulsive model the effective spin coupling 
flows to zero and the spin excitations are gapless. 

The same argument applies to the charge gap at half filling due to the 
appearance of the charge SU{2) symmetry, but the effective couplings are 
reversed; thus a charge gap opens in the repulsive case, but not in the at- 
tractive case. Away from half filling the charge symmetry is broken to U{1) 
and gapless charge excitations are present (holons-antiholons). 

The Z2 particle-hole transformation [1] maps, at half filling, spinous mak- 
ing up the singlet excitation in the repulsive case (always gapless) to the 
objects that make up the the holon-holon excitation in the attractive case 
(gapless at half filling). This is quite remarkable in view of the very different 
structure of the wave functions of the respective states. Not only are the en- 
ergy and momentum the same but also the scattering matrices. To complete 
the corespondance, the triplet of the repulsive model would be matched not 
with the holon-antiholon of the attractive model (which disappears from the 
spectrum at half filling) but with states obtained by removing/ adding two 
electron to the system applying the operator mentioned in Lecure 1. As 
one moves away from half filling one can follow explicitly the breaking of the 
symmetry. 

Spin and charge sectors are coupled; spin excitations involve rearrange- 
ment of charge degrees of freedom and vice versa. In the low energy limit, 
however, we shall see that complete decoupling takes place, similar to the 
one occuring on all energy scales for the Kondo model. 

Let us proceed now to the construction of the eigenstates. We consider 
a system of N electrons moving on a chain of L sites (N and L are assumed 
even, and the lattice spacing is 1 in some unit). We denote by kp = the 
Fermi momentum of the non-interacting system. For the Hubbard model 
aj = sinAjj, and c = |, with m = y, see Lecture 2. Therefore the Bethe- 
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Ansatz equations take the form: 



jk,L _ TT - sm Hj - t 
gj{ A5 - sin % + i 4 

and 

il A, - + if _ -il - sinfe, - if 

}jl Ks-K^-q~ A^ _ sin A;^ + if ' ^ ^ 

for a state with M down spins and iV — M up spins. The quantum number 
M labels both the total spin and its z-component, as discussed earlier, = 
\{N -2M) = S. 

Upon taking the logarithm of eqs(l,2) we obtain: 

M 

Lkj = 27rTCj + ^0(2sin% - 2A5), j = 1, iV (3) 
5=1 

N M 

^0(2A^ - 2sin%) = -27r/^ + ^ 0(A^ - A5), 7 = 1,...,M(4) 

i=i 5=1 

where 



0(a;) = -2tan-^(2a;/'w), -tt < < tt (5) 

The quantum numbers {rij^I^} label the states, and are integers or half 
odd integers. The charge quantum numbers {%} are integers if M is even 
and half odd integers if M is odd. They are defined modulo X, and will 
take values between upper and lower bounds Af^ = — l)/2. The spin 
quantum numbers {/-y} are integers when iV — M — 1 is even, and half odd 
integers when it is odd, and are subject to the same restriction as before, 

r{N,M) = -{N - M -1)12 <I-y< {N -M -1)12 = I+{N,M) (6) 

After solving the equations for an (allowed) configuration {rij^I^} the 
energy and momentum of the that state are 

N 

E = ^-2^cos%, (7) 
i=i 

i=i ^ j 7 
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where the last expression for the momentum is obtained from eq(3)and eq(4). 

The solutions, A-y}, of eqs(l,2) fall, in the thermodynamic limit, into 
patterns described by the string hypothesis [5]: 

1. Real kj 

2. A TC-string, of the form A^f^ = A^") + i^(n + 1 — 2j), z = 1, 2, tc. 

3. A; — A TC-string, consisting of a A ra-string and n complex A;-pairs A;^"^'^, 

each centered around A string position, sin Aj^"^'^ = A^"^ it i^. 

In the Appendix we present a more detailed discussion of some of the 
solutions listed. 

These states form a complete set of 4^ eigenstates [11]. Their degeneracy 
must be properly taken into account. We saw previously that the Bethe- 
states are highest weight states of the spin SU(2) symmetry group, and the 
rest of the multiplet can be obtaines by repeated action of the spin lowering 
operator S~ . Less obvious but still true is the fact that the Bethe states are 
also the highest weight states of the charge SU(2) group and therefore the 
rest of the multiplet can be obtained by the action of the C~ operator. 

We shall be mainly interested in the low lying excitations; the general 
structure of the solutions will interest us only towards the end of the lecture, 
when we wish to sum over all eigenstates to obtain the free energy. 

We turn now to determine which of the possible solutions corresponds to 
the lowest energy eigenvalue. We shall find that the answer depends on the 
sign of the coupling: The ground state is composed of 1-strings (real kj) in 
the repulsive case, and of As 2-strings (rather As — A 1-strings in the terminology 
of the string hypothesis) in the attractive case. Still, the excitations above 
these ground states so differing in structure are spinous and holons carrying 
the same quantum numbers and, in particular at half filling, having the same 
dispersion laws, and interacting via the same S-matrices. 

We consider now each case separately. 

The repulsive Hubbard model. 

The ground state solution consists of real As's and real A's. In the ther- 
modynamic limit we shall be interested in the (real) solution-densities p{k) 
and cr(A), defined by; 

Lp(kj) = !/(%+! 
La{A^) = l/(A^+i - A^) 
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The A;-solutions are distributed between the pseudo charge Fermi- moment a 
— Q and Q with (5 < tt, since they are defined modulo 27r, while the (real) 
A solutions will be distributed between the pseudo spin Fermi momenta —B 
and B < oo. The integration limits are to be determined from 



pik)dk = N/L (9) 

J-Q 

r a{K) = M/L (10) 

J-B 



and can vary with the state, even as the number of electrons is held fixed. 
In terms of the densities the energy and momentum become 



rQ 

E = -2tL / p{k) cos kdk (11) 




P = L p(k)kdk (12) 

J-Q 

The ground state. 

We begin by identifying the ground state configuration /°}. We shall 
find that the ground state is a singlet, Mo = given by a configuration 

symmetrically arranged around zero and as closely packed as possible: 



T° — T° — ^ 
-'7+1 -'7 ~ 

with the rij and filling all the slots from n'^ to n~ and from /"'" to /~ 
respectively, 

I~ < I-y < /+, /± = ±(iV/2 - l)/2. 

n- < rij < n+, = ±(N -l)/2. 

We have tacitly made the choice that N/2 is an odd integer. 

The densities corresponding to this configuration of quantum numbers 
satisfy the integral equations (their derivation is the same as the one outlined 
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in Lecture 3) 

1 fBo 

p^(k) = hcosA;/ dAao(A)Ki(smk - A) (13) 

27r J-Bo 

a,{A) = dkpo{k)Ki{smk- A) - f^" JAV„(A')i^2(A - A') (14) 

J—Qo J— Bo 

Here Po{k), «To(A), Qo and Bo are respectively the densities and integra- 
tion limits for the ground state, and we encounter the kernels from previous 
lecture, 

= - , "\ , = --T^ m , -=1,2... (15) 

with Fourier transform Knip) = e""^!^!. 

The structure of the equations (13,14) will recur often, and it is convenient 
to introduce the notation, 

<^o = ^o + M(Qo,Bo)<^o (16) 

where we denote 

and the matrix of integral operators, 

cos k Ki \ f p 



' i^i -K2 J \a 



(17) 



given in more detail. 



U(n m<b(h W - ( cosk Kii^^^k - A)a{A)dA \ 

^^'^ ^ ' ^ - l^jQ^;^^(gin^_A)p(A;)JA;-/^5i^2(A-A')^(A')iA'j 

To proceed we need to determine the Fermi-levels. We shall find that to 
minimize the energy we have to set Bq = 00 irrespective of the density of elec- 
trons n (which determines the value of Qq.) The state is a singlet as follows 
by intergating eq(14) over A, r^a,{A)dA = ^ = |(f ) = lJ%po{k)dk. 
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Finite values of B correspond to non- vanishing ground state spin and and 
will occur in the presence of a magnetic field. We shall return to this point 
at the end of the lecture. 

As i?o = oo we may solve the spin-equation, eq(14), for any filling, by 
means of the Fourier transform, 

1 fQo 27r 
o-o(A) = - / dkpo(k)sech( — (A -sin A;)). (18) 

U J-Qo U 

and derive an integral equation for po substituting (18) in (13), 

pQ{k) = — + cos A;- Z^" dk'po(k')R f-(sinA; - sin A;')) . (19) 
27r u J-Qo \u ) 

The function i?, the resolvent of i^2, (1 + -R)(l + = 1, is given by [4], 

sech I h'Kt , 



^ ' i-oo 27r 1 + e2|pl 47r J- 



dt- 



+ e^lpl 47r J-00 \^{x-t 



i2' 



It is convenient to introduce the integral operator 
]C(Qo)p(k) = Z^" dk' ls(k - k') - cos k-R f-(sin k - sin k')\ \ p(k'), (21) 

J-Qo I U \U J) 

in terms of which (19) becomes, suppressing the Qq dependence in the oper- 
ator, 

JCpoik) = ^. (22) 

We need to determine Qq. Consider first the ground state at half filling, 
(N = L). The quantum numbers rij assume all allowed values, in other words, 

= Af^ and therefore, also the As-solutions are distributed over the maximal 
range: (^o = ^r, as can be immediately verified in the thermodynamic limit 
by integrating eq(13) over k. In this case we can derive explicit expressions 
for Po{k) and (To(A), 

po(k) = / dp I , 23 

IT Jo 1 + 62 IPI 

(To A = — / dp , ' = / dk sech — A - s n k) 24 

^ ' 27r io ^ cosh(fp) 2Tru .1-^ 'J^ ' 
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Hence the energy can be calculated [1] 



^ = -2t r dkpoik) cos k = -U r j^ y-^^i^) (25) 
L J-TT Jo p(l + 6 2 i^ij 

The ground state away from half-filling is not accessible in closed form, 
except in some limits such as weak coupling, strong coupling or low density 
[4] [12]. In particular, in the strong coupling limit, 1 we have 

4 O 

2Trpo(k) = l + -\n2—cosk, (26) 

U TT 

from which the relation oi n = N/ L and Qo can be determined 

Trn = Qo[l + — ln2smQo], (27) 



and the energy 



Eo/L = -2t 



1 , , 21n2,7rTC,~, sin(27rTC) 
- sin TTTl + — f 1 ^ '- 



(28) 



In the half filled case, = 1, this equation reduces to the ground state 
energy of the Heisenberg antiferromagnet with exchange coupling J = 2t^ jU . 
This relation between couplings can be obtained perturbatively. 

In the general case one needs to solve the equations numerically [4]. The 
results clearly demonstrate the effects of correlations, which become most 
pronounced at strong coupling, or as = N j L tends to 1. 

Elementary excitations 

Let us consider now excitations above the ground state. Excitations are 
created by varying the quantum numbers from their ground state configura- 
tion; spinons - by varying only the spin quantum numbers, {/-y}, and holons 
- by varying only the charge quantum numbers {%}. We shall find that 
holons and spinons are interacting, and will compute their scattering matrix. 
In the low energy limit they decouple in a manner similar to the decoupling 
of holons and spinons in the Kondo model. We shall see that the low lying 
spectrum is linear in the momentum and therefore can be captured by a con- 
tinuum effective hamiltonian: the g-ology model [13] (which is also integrable 
[14].) 
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Spin excitations. 

These are obtained by varying the {/°} sequence from its ground state 
configuartion (the attentive reader will note that the arguments, and words 
parallel those in Lecture 3). 

The triplet. The simplest spin excitation (keeping N fixed) is obtained by 
placing "holes" in the {/-y} sequence. Consider a state with M = N/2 — 1, 
namely a spin-1 state. As before, this corresponds to two "holes" in the 
sequence since the bound I^(N,M) = ±(N — M — l)/2 takes the value 

= ±N/4:. Hence there are N/2 + 1 slots for the spin quantum numbers 
/-y and M = N/2 — 1 of them, leaving two "holes". We shall see again that 
each "hole" corresponds to an excitation, a spinon, carrying spin-|. 

We consider the configuration {/-y} 

leaving the {%} configuration unchanged. (Actually, the rij quantum num- 
bers change from half-odd-integers to integers. This change generates extra 
terms in the energy and momentum which, after a careful analysis, turn out 
not to be relevant.) 

The triplet configuration, eq(29), leads, by methods previously discussed, 
to the following equations for the densities p{k) and «t(A;), 



= if' + M{Q,B)¥ (30) 

where we denote 
with 



a''{A) = S{A-A'l) + S{A-A'i) (31) 

or expilicitly, 

1 /"°° 

p{k) = \-cosk dAa{A)Ki{smk - A) (32) 

2'ir J-oo 

]_ , rQ /-oo 

a{A) = —-a^{A)+ dkp{k)Ki{smk - A) - dA'a{A')K2iA - A/p) 

Ij J — Q J — oo 

We denoted by Aj*, Aj the hole positions corresonding to the quantum num- 
bers omitted from the sequence. To determine them we use the same line of 
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argument as in previous lecture; from the prescribed configurations {rij^I^} 
one finds the corresponding solutions A;j,A-y, and defines the counting func- 
tions. 



M N 



KA) = ^ (E0(A-M-E0(2A-2sin%)) 7 = 1,..., M 

\s=i i=i / 



u,{k) = — [Lk -y2@{2smk -2As)] 
2^ V 5=1 / 



The hole positions in the triplet state then satisfy 

u{Ay = II,. (34) 

Note that the densities p{k) = p{k; Aj^, Aj, Q) and o-(A) = o-(A; Aj^, Aj, Q) de- 
pend also on the hole positions and on the charge pseudo Fermi-momentum, 
Q, which is determined from the condition J^q dkp{k) = ^ and may differ 
from its ground state value Qo- The A integration limit stays at infinity 
since the corrections are of order 1/L. (We shall usually spell out only the 
variables we need.) 

Eqs(32,33) describe a spin-1 state, as can be deduced by integrating with 
respect to A. One finds, f = f + f + ;^(l + l), that is M = N/2 - 1. Again, 
this consideration is just the counterpart in the thermodynamic limit of the 
counting argument presented above. 

Since the holes contribute to the distribution equations terms of order 
we may write p{k) and «t(A) as 

p(k) = poik) + jpi(k) 

a{A) = ao{A) + jar{k) 

with Po{k) = Po{k] Q) and «To(A) = «To(A; Q) satisfying the ground-state 
equations (13,14) for a given value of Q, while pi(k) = pi(k; A'^, A2, Q) and 
«Ti(A) = «t(A; A^, A2, Q) satisfy the equations 

/oo 
JA(7i(A)i^i(sinA; - A) (35) 
-00 

rQ /-oo 
(7i(A) = -a^{A)+ dkpi{k)Ki{smk- A)- JAVi(A')i^2(A -(tSQ) 

J-Q J-00 
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The structure of the equations is as before, $ = + M.{Q, B)^, but now 
the inhomogeneous term is 

^ = ^° + i(-J(A))=^° + i^^ (37) 

and writing $ = $o + -^$1 we derived a reduced equation for $i 

$1 = ^1 + A^(g,5)$i. (38) 

All excitations will be determined from this equation with only the inho- 
mogeneous term Lpi and the integration bounds Q and B varying from case 
to case. 

As before, we solve for by taking the Fourier transform with respect 
to A of eq(45), 

I rQ 27r . 1 gipCA-Af) _^ gip(A-Aj) 

«Ti(A) = — dkpi(k)sech — (A — sin A;) / dp , (39) 

uJ-Q u 2-n J-oo 1 + e 2 IpI 

The second term is identical to the corresponding spin contribution in the 
Kondo model, while the first term represents the rearrangement in the charge 
sector that takes place due to its coupling to the spin sector. 
Feeding eq(39) into eq(35) one has, 

pr{k) = p,iAl A^) = plik- A'l) + plik- A^), (40) 

where each term satisfies, 

;C/9U^;A'') = sech— (sinA;- A'^). (41) 
u u 

We may take the integration limit implicit in the equation to be Qq rather 
than (5, since the correction are of a higher order in Ij L. 

From pi we compute the excitation energy and momentum. The total 
energy eigenvalue is given by 

rQ rQ 
E = -2tL p(k) cos k = Eo(Q) - 2t pAk)cosk, 
J-Q J-Q 

hence the excitation energy, ^E — E — Eq^ 

AEiA'l, A^) = -2t p^{k) cos kdk + [Eo{Q) - Eo{Qo)j , (42) 
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where we define 

fQ 

Eo(Q) = -2tL po(k,Q) cos kdk. (43) 

J-Q 

with po{k^Q) being the solution of the ground-state equation eq(13) with 
integration limit which we now turn to determine. Denoting No{Q) = 
L JSq Poik, Q)dk allows us to rewrite the condition L JSq p{k, Q)dk = N 
as NoiQ) + JSq piik, Q)dk = NoiQo), yielding to order 

Q - Qo = - r Mk)dk(^^^y' . (44) 

•^-Q V OQ y Q=Qo 

Thus to order O(j^) the excitation energy becomes, 

dEo(Qo) 



A£;(A^A^) = -2t [" dkp^(k) cos k + (Q -Qo) „^ 

/Qo rQo 
dkpi{k) COS k — p, / dkpi{k) 
-Qo J-Qo 

where the chemical potential, is given by 

^ = ^= f^MOoh f 9No{Qo) \ -1 ^ _2^_cos Qo + Hk. Qo) cos kdk 



dN V dQo J\ dQo J 1 + J^^^(l){k,Qo)dk 

We denoted 

Mh n\ 1 dpo(k,Q) 

^po(Q,Q) oQ 

and used the fact that Po{k^ Q) = po{—k, Q). 

Again as in Lecture 3 we find that the triplet excitation energy is com- 
posed of two terms /S.E = e(A5*) + e(A2 ) each of which we identify as a spinon 
excitation energy, 

e,(A'^) = - dkpl{k, A'')[2t cos k + p] (47) 

Now to the excitation momentum AP(A5*,A2). The total momentum 
of a state characterized by a configuration {rij^I^} is, as we saw earlier, 

P = ^ Si=i ~^ -^7^ • Here we are considering a configuration with 
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two "holes" in the spin sequence, and hence the momentum of the 
excitation is 

AP(A^ A^) = -^(/f + 4) = -^(KA^^) + KA2)) (48) 

We used Pq = and dropped a term arising upon the shift of all quantum 
numbers rij from half-odd-integers to integers when M changes by one. This 
shift has no effect on the momentum which is periodic in irn. 

We wish to identify the momentum of a single spinon in the sum eq(48). 
It must depend on a single parameter A'* and have a finite limit as X — )• 00. 
The function i/(A'*) implicitly depends on both Aj* and Aj. However, we may 
write 

uiA'l) = J^^ dkQ{2smk-2A'l) - J dA'a{A')Q{A' - A'l) 

= MAl) + ^S{Al,A'^,) (49) 
(and a similar expression for Aj) where, 

y^(A^A^) = (g-g„)^ + y JA;pi(A;;A^A^)0(2sinA;-2Aj) (50) 
L dQo L J-Qo 

-^|jAVi(A';A^A^)0(A'-Aj) 

and Uq is the ground state counting function. Hence, it is clear that the 
spinon momentum can be identified as 

p,{A'') = -—uoiA'') = --n + 2 Z^" JA;/9„(A;)tan-^exp— (A'^ - sinA;),(51) 
L 2 J-Qo u 

yielding, together with the expression for the spinon energy eq(47), a para- 
metric representation of the spinon dispersion. 

The shift of the momentum from a free value, ^(A5*,A2) = L(u(A'^) — 
Uo(A^)), is the scattering phase shift of spinon 1 off spinon 2, according 
to the method of momentum shifts discussed in Lecture 3 [15]. A general 
expression for it is given in eq(50), shall calculate it more explicitly soon. 

We can see immediately that the spinous are gapless; consider e(A'*) in 
the limit A'* — )• ±00. The form of pi(k, A) in this case can be explicitly found 
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from eq(41) [3] 

pUk, K^) = -- cos A;e^^^V(A;), (52) 
u 

with satisfying [4], 

^(A;) = e^'^^"*' + - 1^° dk' cos k'R{-{smk-smk'))ij{k'). (53) 

U J-Q„ U 



Hence, 



p.(A'^) ^ T[|r^-47re±vA''c(°)] 



with 



C(n) ^ ^(n)/ ^ Q^ ^ ]_ f^" ^j.^n.\ cos"(A;) (54) 

27r J-Q„ 



It follows that the energy vanishes linearly with the momentum measured 
with respect to the spin Fermi-momentum kp = = kp, 



e^ip) = ±Vs{p ± k'p) (55) 



with the spin- velocity 



u [c(o) ^ 2^C(o) 



(56) 



These spin-excitations, surviving in the low-enrgy limit, lead to a power law 
behaviour in the spin density correlation function for momentum transfer 
q 2kF. 

The spinous can be excited (in pairs) without exciting charged modes. 
However, due to the coupling of the charge-sector to the spin-sector, the 
charge Fermi sea rearranges when spin modes are excited and modifies their 
energy. This can be seen, in particular, when we study the spinon scattering 
phase shifts. 

To compute the scattering phase shifts we consider the shift of the one 
spinon momentum from its free value ^Z'*, due to the presence of the other 
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spinon, -^/'^ = -^v{K'^) = Ps{A'^) + Thus, scattering phase shift 

for spinons in the triplet state, 

^*"^(A^A^) = -27r(z.(A5^) - z.o(A5^)) 

J-Qo [an J 

/oo 
JAVi(A')0(A'- A^^). (57) 
-oo 

The constant — tt appears because the set {/«} is shifted by | with respect 
to the ground state set. We have 

^*"^'(A5^, A^) = -tt + tIo, 



+ dkpUk-A'l.A'i) |2arctan fe^r^^f-^^"*')) + - + 




with 



The first two terms in (58) correspond to the pure spin scattering and have 
already appeared in the spinon scattering in the Kondo model. The last 
term corresponds to the interaction of the spinon with the change in the 
charge distribution produced by the presence of the triplet. As Qq tends to 
TT, the redistribution of the charge degrees of freedom decreases and so does 
its contribution to the phase shift. 

At half filling the expressions simplify. The Q-level does not shift from 
its ground state value (^o = ^r, so = 0, and the solution of eq(41) is 
straightforward. 



Pi{k) = — cos A; 



— sech — (sin A; — A^) H — sech — (sin A; — At 
u u u u 
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Hence the spinon excitation energy and momentum are 



A'* = — / cosH sech— sinA; - A^)dk = 2t / dp ^ ^ — 

u J -TV u Jo pcosh^p 

pJA'^ = -- + -/ tan-^exp— A'^-smA; =-:t+ / dp ^ . 

2 IT J-TT u 2 Jo pcosh-p 

and the spin velocity at half filling becomes 



2^T7I^, (59) 



with being Bessel functions with imaginary argument. The spinon velocity 
varies from Vg = 2t when = to = when = oo, corresponding to 
the excitation spectrum eg{p) = 2tcosp and eg{p) = 0, respectively. 

The spinon-spinon scattering matrix becomes identical to the Kondo 
spinon scttering matrix (and to the pure spin S-matrix in any SU(2) model 
soluble with the R-matrix given in lecture 2) 

c(tri,u.H .H. _ r(i-f(A^-A^))r(| + f(A^-A^)) 

^'''^ r(i + ff(A5^-A^))r(|-ff(A5^-A^))- ^''^ 



To validate the claim that the spin-1 state considered thus far consists of 
two spin- 1 excitations we must show that a state exists in the spectrum, de- 
generate in energy with the triplet state, with these spin excitations coupled 
to form a spin-0 eigenstate. 

The singlet. The spinous can be coupled antisymmetrically to form a 
singlet by adding to a A- configuration with "holes" at Aj* and Aj a 2-string 
A^ = Ao lb m/4, Ao = KAJ* + Aj). This is, indeed a solution corresponding 
to the choice of /-y-quantum numbers with two unfilled slots at and and 
an additional I^^^ quantum number related to the 2-string position. 

The equations for the (reduced) densities pi,(T\ satisfy the reduced equa- 
tion with inhomogeneous term 



cos k i^2(sin k — Aq) 



(61) 
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where 



a''{A) = S{A-A'I) + ^(A-A^) 
^^*(A) = i^i(A-A„) + i^3(A- A„). 

We find for the spin density 

-Q ^2cosh(fp) l + e-fW 

leading to an equation for pi which is identical to the equation determining 
Pi{k) in the triplet case, (the 2-string contributions, although modifying the 
the spin-sector, cancel in the equation for pi leaving the charge-sector un- 
changed!) Hence the triplet and singlet are degenerate in energy (in the limit 
L — )• oo), confirming the physical picture, totally analogous to the one we 
have found in the Kondo model, of these states describing two spin-| objects 
coupled symmetrically in one case and antisymmetrically in the other. 

The interaction of the spinous, however, depends on the spin-state. Fol- 
lowing the procedure previously outlined we find that the spinous, when in 
the singlet state, undergo scattering with the phase shift 

S^ing^^h^^h^ ^ (62) 





^f(A5^-A^)J i " -AL^^ra. -AizAi 



2 



+ 



r dkp,ik;AtA',) |arctan(ev(Af-^in.)^ _ - _ 
J-Qo [ 2 dn ) 



consisting of a charge contribution identical to the one we encountered dis- 
cussing the triplet scattering, and a spin contribution that is modified in a 
manner similar to the modification we found for the Kondo-spinons scattering 
in the singlet state. 

Again, in the case of half filling the charge sector does not contribute to 
the spinon scattering, and the singlet S-matrix is found from the first two 
terms on the right hand side of eq(62) 

.isin^^^.H _ 1 + m - A2) r(i - ^(A^ - A^))r(| + ^(A^ - A^)) 
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We can combine (60) and (63) to get the scattering matrix for spinons at 
half fiUing 

_ r(l - ^(A^ - A^))r(| + ^(A^ - A^)) f (A^ - A^)/^^ + zfP^^ ] 
^ " r(l + i(A5^-A^))r(|-i(A5^-A^)) 1 (A^,-Al) + q J 

where I^^ and P^^ are the identity and the exchange operator in spin space, 
respectively. This expression satisfies the Yang-Baxter equation guaranteeing 
that the physical, dressed S-matrices which describe the scattering of the 
actual quasi-particles do factorize consistently in the same manner as do the 
bare S-matrices discussed in Lecture 1. As a consequence we conclude that 
the excitations, though interacting, never decay! This is more remarkable still 
since gapless excitations are present, and is due to the dynamical conservation 
laws (briefly mentioned in Lecture 1) which protect the excitations. 

Charge excitations. 

Still keeping the number of electrons fixed, we now vary the {%} se- 
quence from its vacuum configuration, leaving the spin configuration {/-y} 
unchanged. 

The holon-antiholon. We consider first the case N < where the fol- 
lowing charge configuration is allowed; remove the level ha from the ground 
state sequence (creating a hole at the corresponding A;^ < Qo), and add a 
level outside the charge-sea, hb (creating a particle at kB > Qo)- Obviously, 
this excitation is not present at half filling. 

The {%•} configuration we consider, hb > Ji"*", %+i ~ = 1 + ^j,A^ 
leads to the equations (here —Q<k<Q) 

11 /"°° 

p[k) = —S{k — kA) + cos k dAa{A)Ki{sm k — A) 

2'ir L J-oo 

1 rQ /-oo 

a{A) = -Ki{smkB-A)+ dkp{k)Ki{smk - A) - JAV(A')i^2(A - A 

IJ J —Q J — oo 

with the Q-level set by the requirement JSq p(k)dk = It can be easily 

checked integrating over A that M = |iV as expected since we did not change 
the spin-sequence {/°}. 

Again, pi{k) and cri(A) satisfy the reduced equation with the inhomoge- 
neous term 
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Hence immediately, 



fQ 

^i{p) = / dkpi{k) 
J-Q 



^—ipsink ^—ip sin kj^ 



and the equation for pi{k) follows. It is convenient to introduce a smooth 
density p'i{k)^ 

p[{k) = p^{k) + S{k-kA) (66) 

and one has 

p[{k) = p[{k- kA, ks) = plik- kA) - pl{k- ks). (67) 

with 

4 /4 . . \ . 

]Cp\{k]kj) = cos A; i?( — (sin A; — sin A;j) j j = A,B. (68) 

The calculation of the energy and momentum proceeds as before; the 
energy is given by 

fQ 

E = —2t cos kB — 2tL / dkp(k) cos k 

J-Q 

= -2tcoskB-2t dkip'Ak) - 8ik- k a)) cos k + EoiQ) 

J-Q 

and the excitation energy by 

rQ 

AE^kA^kB) = 2^ cos A;^ — 2^ cos A;^ — 2^ / dkp[{k; kA^kB) cos k -\- Eo{Q) — Eg 

J-Q 

= -^c(kA) + ^c(kB), 

where we introduced the holon energy 
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ec{kA) = —2t cos + / dkpKk^ kA)[2t cos k -\- fj]. (70) 

This function is monotonically decreasing to 

ec(Qo) = (71) 

as kA tends to Qq. 

The momentum of the excitation is 

27r 27r 27r 

P - Po = — XK^' ~ + 'Y^^ ^ ^^^^ ~ ^^"^ ^ -i'cl^^) + Pc{kB), (72) 
i 

with the holon momentum, following previous arguments, 

VcikA) = -T'^oikA) = 27r / po{k)dk (73) 
L Jo 



i: 



dkpoik)log ' ^ « y V2 



sinAj/i — sinA: \ "p / 1 ' sink a —sink 



p / ^ sinfc4 —sin k \ f 1 i ^ sinfc4 — sinfc 

* u / \ 2 u 



Both forms for the momentum will be useful. Together with the expres- 
sion for the holon energy we have a parametric representation of the holon 
dispersion. 

It is obvious that the excitation is gapless when kA = ks = Qo (recall that 
— TT < —Qo<kA 1^ Qo < < TT.) The holon momentum in this limit tends 
to the charge Fermi momentum kp = irn, Pc(kA) = 27r J^^ po(k)dk — )• zLtttc, 
hence the vanishing of the excitation energy will lead to power law behaviour 
for charge density correlation functions at momentum transfer q 2kp = 
Akp- In a few paragraphs we shall find that the holon-holon excitation, on 
the other hand, is gapful. 

Once more we are able to describe the excitation as a combination of two 
objects. We will call holon the one with energy ec(A;^) and momentum Pc{kA) 
associated with the hole in the sequence, and antiholon the one created by 
adding an electron to the system and characterized by —ec{kB) and —pc{kB)- 

The holon- antiholon phase shift = —2'K(u)(kB) — <^o(^b)) is given by 

1 f r fl + ^^iiLiaii^iiL^^ Y (- — ^^iiLiaii^isJA 

J- ^ I \ ' 11. 1x9. 11. 



^ ^ I r I 1 ^ sinfcR— sinfc4 j p | i _|_ ^ sin r — sin / 
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/ dkp[{k; 
J —Q(\ 



As before there are two contributions to the phase shift: the first Hne cor- 
responds to direct holon-antiholon interaction; the rest describes the effect 
of the interaction between the excitations and the redistribution of charge 
they produce. Since the interaction is between charge degrees of freedom, 
the integrand is of the same form as the first term of the phase shift. At 
half fiUing it vanishes and the first term would describe the holon-antiholon 
scattering had they existed. 

The holon-holon excitation. We were considering configurations with only 
real A;-momenta and holes in them. To discuss states with double occupancy, 
however, we need to consider solutions involving complex A;-pairs. As a result, 
a gap of order U opens, and the excitations lie in the "upper Hubbard band". 
These gapful excitations are present both at half filling and away from it, but 
are the only ones surviving at half filling ( keeping the number of electrons 
fixed.) The disappearance of the gapless holon-antihilon excitation from the 
spectrum at half filling is the origin of the charge gap and the concommitant 
metal-insulator transition. 

Consider then a configuration {%} leading to a complex A;— pair k^ = 
K i ix and two "holes" in the real A;-sea (the number of electrons is kept 
fixed). In the notation of the Appendix we have M' = 1 complex A;-pairs (i.e. 
a A; — A 1-string) with the associated quantum number /'. The number of 
real A;-momenta is iV — 2 and M = N/2 — 1. The number of real A decreased 
by one since one of them becomes the parameter A to describe the complex 
pair. From the counting arguments we have /' = 0, the set {/-y} splits around 
zero so that J^-y -^7 = 0, but there is a tt contribution to the phase shift. 

The Bethe Ansatz equations now take the form 



M-l 

L{k ± ix) = 2Trnj + ^ 0[2 sin(/€ ± ix) - 2As] (75) 

5=1 

The right-hand side can be converted to an integral and evaluated to 
order by means of the ground state density (To(A). One finds the complex 
is driven to string positions centered around a particular A solution. Denoting 
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this particular solution by A, we find that it must satisfy (to order 0(e 
with 1] a number of order 1) 

sin(/€ ± ix) = A T W4 + 0(6''^) (76) 

namely, 

u 

k± = TT — arcsin(A it i—) (77) 

with A real. 

Substituting (77) in the eigenvalue equations for real A;, 
^ik,L ^ A^-sin%-if \ A - sin % - if 

it follows that the real A;'s (the 1-strings) satisfy 

M-l 

Lkj = 2Trnj + Y, ©(2 sin % - 2As) + 0(2 sin % - 2A) (78) 

5=1 

where the {%} set has two "holes" in it, at n'^ and TCj, corresponding to 
omitted momenta and , respectively. 

Similarly, the eigenvalue equation for real-A, becomes 

_ 1^-' As-A-y + i^ \ A - + if _ (^TT^ A-,-sin%-if \ A-, - A - if 
~ A5 - A^ - if j A - A^ - if " [fJi A^ - sin% + if J A^ - A + if 

namely, 

N-2 M-l 

^ 0(2A^-2sin%) = -27r/^+ ^ 0(A^- A5), 7 = 1,...,M-1 (79) 

i=i 5=1 

which has the same form as eq(4) but for N — 2 particles. 

Solving the equations for real k and A and evaluating eq(75) one has 

A = ^(sinA;f + sinA;^) (80) 
similar to the result we encountered discussing singlet spin excitations. 
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We convert (78) and (79), to the integral equations for p{k) and cr(A), 
then defining pi and cti, and finally introducing 

p[{k) = p,{k) + jS{k-k'l) + jS{k-k'^) (81) 

we find that the densities ai and p[ satsfy the usual reduced equation with 
the following inhomogeneous term 

_/ cos A; i^i(sinA; - A) \ , . 

" \-K^{smk'^ - A) - K^isink^ - A) J ^ ' 

and with the Q-level set by the condition J^q p{k)dk = ^^yT' 
Hence, 

Mp) = I_^pW2,oMIv) ^^^) ^^^^ 

and 

p',{k) = p[(k; kl k'i) = p^(k; k'l) + p^(k; k'^) + p'\k; A) (84) 
where p''(k;kj) satisfy eq(68), and /9**(A;;A) is the solution of 

A:/9'*(A;;A) = cosA; i^i(sinA; - A). (85) 
The energy of the state is given by 

E = -2t{cosk+ + cosk_) -2tL dkp{k) cos k (86) 

J — Q 

with the first term being the contribition of the string. Using some identities: 



cos A;+ + cos A;_ = 2Re [^1 - {X + i^Y) = ^ + ^ J^^ dk cos^ kKi{smk - A), 
we find for the excitation energy 

AE(k'l, k^) = U + J2 ^c(k^) + 2t r dk cos^ kK^ (sin A; - A) 

i=i 

-2t Z^" JA;/9f(A;;A)|cosA; + -^| 

= A,(U,Qo) + j2e,(k';) (87) 
i=i 
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with the ec{k) same as defined for the holon-antiholon case, eq(70). The 
holon-holon energy gap Ag{U^Qo) + 2ec{Qo) does not vanish for any fiUing. 
Away from half fiUing, however, the holon-antiholon gapless modes are avail- 
able to carry the charge. 

We evaluate now the excitation momentum. Since there are no holes in 
the spin quantum numbers and we are taking = 0, the only contribution 
to the change in momentum comes from the holes in the charge distribution. 
Then 

AP = -^{nU4) (88) 

f -1 

= k'l + k^-j J2 {0(2(sin k'l - Ap j) + 0(2(sin k^ - 

^ 13=1 

+ y {0(2(sin k'l - A)) + 0(2(sin k^ - \ 

= pM)+Pc{k^2)- (89) 

with Pc{k) defined in (73). 

The holon-holon phase shift is 

\2p„(Q„) + /_* dQ„ 

/OO 
dAai{A)@{smki - A) - 0(2(sinA;i - A)), 
-OO 

and can be written more explicitly as 
,,,,,, 1 l + i^(sinA;'^-sinA;,^)r(l-^^^^ ^)rU + ^ ^ ^ 



t I 1 — (sin ki — sin ) p _|_ ^2 ^ p | 1 



+ 



dkp[{k;klk'^) 

'J —Q(\ 



t T I -I -siak^-siak \ -p / 1 , -siak^ -siak \ (in 

r (^1 - ^ — — ) ' — — j 



It is worth comparing (74) and (91), specially the terms that do not depend 
on the charge densities. These, as we shall see, are the only terms surviving 
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at half filling. If both excitations were present at half filling, they would have 
the singlet-triplet relation that we found for the spin excitations. 

At half filling matters simplify again; we have explicit solutions 
4 4 

p^ik; kf) = — cos k-R(-(sm k — sin kf)), p'^Hk; A) = cos A;i^i(sin k — A), 
u u 

leading to the holon-holon excitation energy 

I^E'^^^k'l.k'i) = U ^j^e^ik';), (92) 

i=i 

with the holon energy explicitly given by 

e:ik'^) = -2t cos k'^ + At (93) 
and the gap by 

A, = [/ + 2e:(7r). (94) 
Likewise, the momentum of the holon can be explicitly computed 

I'k'' /-oo p-fP Jn(n) 

pAk^) = 27r / po{k)dk = k^+ dp -^sin(psinA;'^) 

Jo Jo p cosh ^p 

p I _j_ ^ sinfc^ — sinfc^ \ ^ sin fc^ — sin k' 



1 r M- J- 

k'' + : / iA;'log <^ — ^ . . . ,x / 

27rZ J —ir I r I 1 ^' sin fc" — sin k \ ( 1 \ ^ sinfc" — sm 



and the phase shift, 

{-n I -I -sinfco — sinfcj' \ -n / 1 i • sinfc? — sin fc? \ 
1 - i^(sin A;,^ - sin A;,^) ^ ( 1 " J ^ U + ) 
-2 • A ^^—7 ( 7 
1 + t-[sm — Sm ki) p | ]^ _|_ ^ sinfc^-sinfcf J p | 1 _ ^ sin fc^-sin fcf 

Note the similarity with the singlet scattering shift, eq(63). 
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Charge-spin excitations. 

We studied thus far pure spin excitations, as well as pure charge excita- 
tions. The latter fall into two categories; the gapless holon-antiholon and the 
gapful holon-holon excitations. In the half filled band only the latter exists 
and a charge gap opens. Similarly the simplest spin excitation are the triplet 
and the singlet composed of two spin-1/2 spinous. 

A single holon cannot be excited if we keep the the number of electrons 
fixed. Neither can we excite a single spinon. To do so we need to change the 
number of particles in the system. When we add an electron we modify both 
the spin and charge configuration, creating a hole in the former and adding a 
level in the latter, thus exciting a spinon and an antiholon excitation (rather 
a coherent superposition of them). When we remove an electron we create 
holes both in the spin and the charge sequences, exciting a spinon and a 
holon. These considerations are similar to the corresponding ones in the 
Kondo model. However, unlike the situation there, the spinon and holon do 
not decouple except in the low-energy limit. 

Now that we are considering a change in the number of particles it is 
convenient to change the form of the hamiltonian to make it more symmetric. 
We shall replace the interaction term Hj = U Y^j by Hj = + 

Tiji — ly. The difference, H'j — Hj = —^UN = —\U + commutes 

with the hamiltonian and only modifies the energetics. We have, in other 
words, added a chemical potential, A = —^U, to the system to make it 
particle-hole symmetric. 

We begin by removing an electron, creating a charge hole at k'^ and a 
spin hole at A'* corresponding to the omitted quantum numbers tc'*, J'*. The 
rij quantum numbers change from half-odd-integers to integers, while the 
stay half-odd-integers. 

Following the well trodden path we introduce the densities pi{k) and 
«Ti(A), then p[{k) = pi{k) + S{k — k^) obeying 

1 27r 4/4 
]Cp[(k) = — cos A;— sech( — (sin k — A'*)) — cos k—R ( —(sin k — sin k'^) 
u u u \u 

with the solution given by 

p\{k)=p{{k-K^) + p\{k-k^). 
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Hence the excitation energy 

E(N - 1; A'^) - Eo(N) = AE^k'^, ^''^ = J + ^^(^'^^ + (^6) 

and momentum 

AP = (/'^ + n'^) = p,{k^) + p,{K^), (97) 

Therefore, this excitation is composed of a holon and a spinon. They interact 
and the phase shift can be calculated, for instance, from the spin contribution 

to be, 

ch,s ^ i ( — (A''-sinfc'')\ 

' = z arctan e ^ ^ ' ) 

- l"^" dkp[(k) |2arctan(ev(A''--fc)) _ _ ^1 ^gg^ 
J-Qo [ 2 dn ) 

When we remove an electron in a definite momentum state p forming the 
state Cp^al^ > we end up with a coherent superposition of the eigenstates just 
described, subject to p = pdk'^) +Ps(A'*). The associated spread in energies 
gives a measure of the life time of an electron (or electron hole) in the system. 

At half-filling the excitation energy and momentum are 

. ^ /,?, .?,x U , ^ Ji(p) cosipsink'^) , Ji(»)cos»A'' 

AE.ik^A'^) = - + 2tcosk'^ + U dp ^^-^ ' +2t dp '^"^^ ^ 

2 Jo p[l + e2'P') Jo pcosh-p 

AP{k\K-) = r^/o(p)sm(M'^) 
2 Jo 



p cosh ^p 
+k^ + / JA;'log<^— ^ 



4 

siuk^ —siuk' \ "P / 1 * sin A:'* — sin A:' 



r (1 + -smfe ) r (i 

: sin —sink' \ Y" ( — -L -i sin fc^ — sin k' 
^ I 2 * u 



Jo pcosh-p Jo pcoshp- 



The minimum energy to remove a particle from a half filled band occurs 
at k'^ = TT, A'* = cxD, 

1 /i (v) 1 1 

AE_ = E^{N-l)-E^{N) = -U-2t+4t / dp '^Z.. = ::U+e,{7r) = -A,(7r). 

2 Jo » ( 1 + ealPl ) 2 2 



p ( 1 + 62 I 

(99) 
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We shall compare it to the energy required to add a particle at half filling in 
the next subsection. 

The phase shift is just the first term in (98) 

Sf^,s ^ ^ _ 2arctan(ev(A''--fc'')^^ (100) 

and is similar to the spinon-impurity scattering phase shift we encounterd in 
the Kondo model. The reason for the similarity is that in the half filled case 
in the Hubbard model the charge distribution is locked and does not interact 
directly, while in the Kondo model it decouples completely. We conclude 
that spinon-holon scattering matrix 

5;--(fc\A-) = i ,+ (101) 

indicates that the low-energy spinon sector (A — )• iboo) completely decouples 
from the charge sector. This decoupling is captured by the g-ology model, 
describing the low-energy physics of the model. As discussed earlier, in the 
latter model spinon sector and holon sector belong to different Hilbert spaces 
in analogy to the situation in the Lecture 3. It was in this context that the 
first spinon S-matrix was calculated [15]. 

Consider now adding a particle to the system. This can be done in two 
ways; creating a state with double occupancy, described in the Bethe-Ansatz 
language by configuration containing complex A;-pairs, or creating a state 
with all momenta real. At half filling only the first possibility exists. 

We begin by discussing this situation. We consider a configuration with 
one hole in the real A;-sea and a 2-string. A hole also opens in the A-sequence 
since N is increased by one while M is unchanged. Placing the charge-hole 
at A;'*, the charge-string at A;^, where sin A;^ = A it z^, and the spin-hole at 
A'*, we find upon solving for the density of real-A and real-As that A = sin As'*. 
The relevant integral equation is 

f 1 27r 4/4 
KpAk) = cos k \ i^i(sin As — A) sech( — (sin As — A )) R I —(sin k — sin 

leading to the excitation energy, momentum and phase-shift 
AE^^\k^,A'') = E^^\N + l;k^,A'') - EoiN) 
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= A,([/,go) + e.(A'^) + ec(A;'^), (102) 
AP(k^,A^) = -^(lf^ + n^)=p^(k^)+p^(A^), (103) 
Ss,m^f,h^^h^ = --2arctan(e^(^''-^^"*"')) 

- l"^" dkp[(k) |2arctan(ev(A.--fc'')) _1_ ^1(104) 
J-Qo [ 2 dn f 

At half filling we have explicit results, and it can be easily seen that they 
are the same as when we were removing an electron, 

l^r , rh , r J Ji(p)cospA'' , p , Ji(p)cos(psinA; 



AE^\k^, A^) = ^U + 2t cos k^ + 2t H + At H dp'^ 

2 Jo p cosh jp Jo 7,(1 



p[l + e 

where is composed of two contributions: \U — the chemical potential, 
and U — string contribution arising from double occupancy. Hence the min- 
imum energy required to add an electron, 

E^{N + l)-Eo{N)='^U + e,{T,) (105) 

is the same as to remove one (due to our choice of a particle-hole symmetric 
hamiltonian). 

Away from half filling we can add a particle without creating double 
occupancy. There are spaces available for an extra k = k'^, but as before a 
hole opens in the {I-y} sequence at position A'*. 

The excitation energy is 



AE^^\k^,A'') = -2t cos k^ - dkpi{k){2t cos k- n} 



u 

2"' 



-| + e,(A'^) + e,(A;'^), (106) 



and the momentum 



AP(A;^ A'^) = ^in^- l"-) = Pc{k^) + p.(A'^). (107) 



L 

The spinon-antiholon phase- shift 



- + 2arctan(ev(A''--fc'')) _ |arctan(ev(A''--fc)) _ _ ^ 

2 J-Qn 2 dn 
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is conjugate to the spinon-holon phase-shift. 

To summarize our conclusions for the repulsive interaction: the spin ex- 
citations are gapless for any value of tc, while gapless charge excitations exist 
only away from half filling. When the density reaches this critical point, all 
gapless charge carrying modes disappear from the spectrum (gapful modes 
are always present), and the system becomes an insulator. This is the mech- 
anism underlying the Mott transition in this model. 

The structure of the excitations is rather complex. Away from half filling, 
there is an incoherent charge background that is modified in the excited 
states. This is particularly manifest in the phase shifts. We see that, for 
a given excited state, the phase shift of one of the elemantary excitations 
consists of a term due to the interaction with the other excitation, and a 
term that describes the interaction with the change from the ground state of 
the charge incoherent background. 

At half filling, the charge distribution becomes rigid^ and we only see 
interactions between the elementary excitations. 

When a magnetic field is turned on similar comments will apply to a spin 
background. We will see the same behavior in the attractive case. 

The low energy physics can be captured by an effective fixed point hamil- 
tonian, calculated by repeatedly integrating out higher energy degrees of 
freedom. This effective hamiltonian in our case will describe excitations with 
linear dispersion and hence will be conformally invariant. One can, then, 
from the exact solution calculate the parameters specifying the conformal 
hamiltonian directly without carrying out RG transformations [20] [21]. 

The attractive Hubbard model 

The structure of the wave functions of the ground state and the excitations 
changes completely; it is advantageous to have doubly occupied sites to lower 
the energy. Therefore the ground state is composed entirely of k 2-strings, 
and the excitations involve either the breaking or displacing of 2-strings. 

The ground state 

The ground state configuration for u < consist of a sea of k 2-strings. 
To minimize the energy we choose the configuration of quantum numbers 
to consists of consecutive integers centered around 0, filling all slots 
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between / and / . The k 2-strings satisfy the string hypothesis 

u 

ks,± = a.T:csm{As ± i-) S=1,---,M, (108) 

where Ag are real solutions of the eigenvalue equations. It is easy to see that 
this choice satisfies the A-equations (2) trivially, whereas the A;-equations 
become 

/Ma_a_-u\ 

whose logarithmic version is 

\u\ ^ 
2i:i?e{arcsin(A„ + )} = 27r/:, - ^ 0(A„ - A^j)), (110) 

4 13=1 

with integer (h.o.i) if X — M is odd (even). The integral equation satisfied 
by the density of strings, (To(A), for the ground state in the thermodynamic 
limit is 



jCaoiA) = aoiA) + f^" JAVo(A')i^2(A - A') = -Re \ . \ 

J-Bo TT ( Jl - (A 



■111) 



where now Kn(x) = - — , and the integration limit Bq depends on the 

density n = N/ L through /_^^ dAao(A) = ^N/ L. The integral operator C 
introduced in eq(lll) plays a role similar to the one played by the operator 
}C in the repulsive case. 

The energy and the momentum of the ground state are given by 



" /'-DO / 

Eo{Bo) = 2tJ2i^os ks,+ + cos kg,-) = -UL dAa{A)Re\ 1 - {A - 

5=1 J-Bo V 4 

N „ /--Bo 

= -\U\ 2t dkcos^k dAao(A)Ki(smk - A), (112) 

2 J-TT J-Ba 



9 M' 

a = l 
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where we explicitly separated the contribution to the energy, — -y^iV, due to 
double occupancy. 

Consider first the case of half filling. It is obtained by setting i?o = oo, 
as can be easily verified by integrating eq(lll). In this case the equation can 
be solved by Fourier transform to yield 

^^^^^ = ^1 cosh (M,^ ' ^'''^ 



the same function that described the ground state spin density in the U > 
case, eq(24). From eq(113) and eq(112) we get 



JoMAM_. ,114) 



Comparing (114) to the results for the U > case, we find a relation between 
the the ground state energies of the attractive and repulsive hamiltonias 

<*([/; 5, = oo) = ^ + E^,^^i\U\- Qo = tt) 

as expected on general grounds [1]. 

To go away from half- filling one needs i?o < oo, and it is convenient 
to use a formalism where the physical quantities at arbitrary density are 
given by equations with those quantities evaluated at half filling appearing 
as inhomogenous terms. The method is due to Griffith [16]: eq(lll) is of the 
form 

a{A) = /(A) + Z^" dK'a{K')K{K - A'), (115) 



and from the solution we wish to calculate integrals of the form 



dKa{K)eo{K), (116) 

J-Bn 



where /(A), K{K — A'), and eo(A) are known functions defined in the whole 
real axis, with Fourier transforms /(p), K{p), ^{p)- this particular case 
/(A) = </.o(A) = ii?e|l/^l-(A + iM)2| and eo(A) = Re^l - (A 
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The function «t(A) is physically relevant in the interval [—Bq^ Bo]. How- 
ever, since /(A) and K{A) are defined in the whole real axis, the integral 
equation (115) defines a continuation of «t(A) to the whole real axis. There- 
fore, we may manipulate the Fourier transform in (115), 

Hp) = f(p) + r d^'H^Y'^'Kiv). 

J-Bo 

[J-oo J\A'\>Bo ) 

= Hp) + K{p)Hp) - I dA'a{A')e'^^'K{p) 

^|A'|>Bo 

to obtain an equivalent equation for cr{A) 

(t(A) = (t^(A)- f dA'(T(A')R(A- A'),. (117) 

^|A'|>B<, 

The function 

^-(A) = — r dp^^^^S^ (118) 



is the solution of (115) for i?o — ^ oo, and we have often encountered the 

resolvent R(A - A') = ^ dP^^'^^l^- 
We rewrite now the integral (116) 

dAa{A)eoiA) = I - \ dAa{A)eoiA) 

-Bo [J-oo J\A\>Bo) 

/CO 
JA(7°°(A)eo(A) 
-QO 

- / / dAa{A')R{A - A')eo(A) - / dAa{A)eoiA) 

J\A\>Bo J-oo J\A\>Bo 

= /°° - / dA(T(A)e(A), 

J\A\>B 

where e(A) is given by 

/oo 
dA'R{A - A')eo(A'), (119) 
-QO 
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and /°° is the value of / when Bq ^ oo. 
In our particular case we have therefore, 



J^ao{A) = ao{A) - A / dA'ao{A')R 



\u\ J\A\>Bo 




<W (120) 



where is the density at half filling and is given below. The integral 
operator J- will recur in our subsequent discussion. (This type of equation 
also describes the magnetization of the repulsive Hubbard model, with Bq 
related to the magnetic field.) 
The ground state energy is 



EoiBo)/L = E^/L + \U\SiBo)/2 + / JA^o(A')er (A'), (121) 



where S = S{Bo) = l — n{Bo) measures the doping, and is the the function 



we met earlier when discussing spinon excitations of the repulsive model at 
half filling, = ej. We shall see soon, when we discuss charge excitation 
in the attractive case, that it also corresponds to the energy associated with 
the charge excitations, hence the notation. The fact that the same function 
appears in the ground state energy and in the excitations is due to the struc- 
ture of the ground state in the attractive case. Eq(121) allows us to view the 
ground state energy as the energy at half filling plus the energy of the holes 
that have been made on the charge distribution in order to reach a particular 
filling. 

Close to half filling we find an expansion in ^, 



|A|>Bo 






(123) 



From this, the chemical potential is obtained 




(124) 
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(Notice that the chemical potential remains negative for any \U\ ^ 0.) Fi- 
nally, we write the charge susceptibility 

dn 1 -^0 (o 



^_aHa. (125) 



Elementary excitations 

To obtain excitations one has to consider small variations of the quantum 
numbers away from the ground state configuration. This can be done with 
pair-breaking (spin excitations) or pair-rearrangement (charge excitations). 
We shall discuss both possibilities. 

Charge excitations 

The Holon-holon (Quartet). This excitation exists at half-filling and away 
from it. It is gapless and resembles the spinon excitation of repulsive model. 
To create it remove two 2-strings and combine the resulting four momenta 
into the following configuration, 

I I 

kl = arcsm(A:: + z-), = . - arcsm(A::), 

I ^ I 

= arcsm(A::), = arcsm(A:: - z-). (126) 

The quartet is parameterized by A" and the quantum number associated is 
/". 

a 

From the constraints on the quantum numbers we can see that, at half- 
filling, /" = 0, and the set is not shifted with respect to the ground 
state configuration. We will assume that this holds away from half-filling. 

The set A^ is determined therefore from, 

2i:i?e{arcsin(A^ + i^)} 

N/2-1 2 \ 

= 2<- E 0(A:,-A;,)-|0(2(A:,-A")) + 0(-(A:,-A"))|i27) 

which translates in the thermodynamic limit to an integral equation for the 
k 2-string distribution, a{A;B)^ 

Ca(A) = </.o(A) - I (^(A - A^^) + ^(A - A^) + K^(A - A") + Ks(A - A")) (128) 
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with A" = |(A5* + A2). The integration Hmits are set by the condition 
LJ^gdAaiA) = N/2 - 1. 
Writing a{A;B) as 

^(A; B) = ao{A; B) + ^a[{A- B,, A^ A^) - ^{8{A - A'l) + ^(A - A^)) (129) 

we find that the excitation density is given by 

^;(A) = ^^(A;A5^) + ^^(A;A^)-^?(A) , (130) 

where 

Cal{A-A^) = i^2(A-Aj), (131) 
C(tI(A; a'') = i^i(A - A") + Ks(A - A") . (132) 

These equations bear similarity to the singlet equations in the repulsive 
case. Here, as note before, they are defined with respect to finite integration 
limits ibi?o while there, as long as we considered excitations out of the ground 
state, Bq = 00. Rewriting the equations as, 

^^^(A); A'^) = (A; A"), ^^?(A; A^ A^) = a^A; A"). (133) 
we find 

^^quartet^^h^ A^) = e,(A^) + ^A'^) - eg(A") . (134) 

with 

MAJ) = C(AJ)+ /, „ '^A^i=(A;Aj)(M + ^ + ,oo(A)\ 

J\A\>Bo [ ^ ) 

e,(A") = / JA^7?(A;A"){|[/| + 2// + C(A)}, 

where e^(A'*) the charge excitation energy at half filling, to be evaluated 
shortly. 

The chemical potential fj, is 

dEo y^-iBo + i^-^y] + I\ iA*o(A)i?e{^l-(A) + iM)2| 

dN 1 + JX JA*o(A) ^ ^ 
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with *o(A) = satisfying /:{*o(A)} = -^{^^(A - Bo) + MA + 

Bq)}. Combining this result with the expression for one finds = 
Ag(i?o) is a gap similar to the one encountered in the repulsive case, with the 
important difference that it vanishes at half filling, Ag(oo) = 0. We shall see 
next that holon-antiholon excitations are gapless but, of course, disappear 
from the spectrum at half filling. This is precisely where the holon-holon 
ecitation becomes gapless. We conclude therefore that in the attractive case 
gapless charge excitations are always present. 
The excitation momentum is 

^p,uartet^^k^ A^) = "^(^(Al) + H^2)) = Pd^'l) + Pc(A^), (136) 

with the counting function u defined as usual (in this case from eq(127)). 
The momentum function pc(A'*) = —^Uo(A^) takes the form 

p^{A'') = -2i?e{arcsin(A'^ + i—f} - 27r / JAVo(A')0(A'^ - A') 

4 J-Bo 

= p^{K^) - i / dk'aUK'- K^) log \ —) "I ( y . 

In terms of these functions the excitation energy and momentum of the 
quartet state are 

AE'^uartet^^h^^h-^ = e,(A5^) + e,(A^) + A,(5o) 
Api^-rtet^^h^^h^ = p,(A^) + p^A';) (137) 

As before, the quartet can be described as the combination of two objects, 
in this case both are holons. They interact and their scattering is given by 



^quartet ^j^h J^^^ 



1, |l + ^^(A5^-A^)^^l-^-^j^^i + ^ 



Tlog 

I 
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where 



. #.(Am) _ Q(A^ - Bo) + 0(A^ + Bo) + J\ dA'^o(A')Q(Al - AQ 

dn l + /\JA*o(A) ^' ^ 



We now turn to compute the quantities at half filling. As i?o — ^ oo we 
have explicitly, 

(7^'°°(A) = ^i? ( ^(A - A^) ) , and (7?'°°(A) = i^i(A- A"). (139) 
\u\ / 

Hence 

e;f (Aj) = 2t r dkcos^kKi{smk- A) 

J —IT 

''0 T) cosh — T) 



pcosh y^p 

whereas = 0. The function e^(A) was arrived at by a different route 
when we discussed the ground state, eq(122). 
We also have 

and the chemical potential becomes 

f^ = -\U\/2. 

We observe that at half filling the holons of the attractive model have 
the same dispersion as the spinous of the repulsive model. This is quite 
remarkable in view of the profound differnce in structure of the respective 
ground states. It is, of course, due to the underlying SU{2) X SU{2) symmetry 
that is manifest at half filling, and the Z2 transformation that exchanges 
them. 
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The phase shift at half filling follows immediately, 

(142) 

Comparing this result to the spin excitations for the U > case, we note 
again, 

substantiating the identification of holons in the attractive model with spinous 
in the repulsive model. 

Holon-antiholon. This excitation is only present away from half-filling, 
where it is possible to remove one pair with A'* < Bq from the sequence, and 
add another with A'* > i?o, namely outside the {A^} set of the ground state. 
Denoting by j'*, j'* the corresponding quantum numbers, the excitation mo- 
mentum becomes 

^ph,h ^ ^^jh _ jh^^ ^^44^ 

and the integral equations for the A-density 

Ca{K) = </.o(A) - |^(A - A'^) + ^K^^K - A'^). (145) 

Again, writing cr(A) = «To(A') + ■^«t^(A) — |;^(A — A'*), the integral equation 
becomes 

/:(7;(A) = i^2(A- A'^) -i^2(A- A'^). (146) 

with the solution 

a[{K) = (tI(A - A'') - (tI(A - A^), (147) 

leading to the excitation energy and the momentum 

l^E^,h{A^,A^) = e,(A'') - e,(A^) 

l^P^^^{A^,A^) = p,{A^)-p,{A^). (148) 



^log 

I 



r 1 



-A'' 



A^) 



Aft 
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The holon-antiholon excitation is therefore gapless at A'* = A'* = Bg. At 
this point the holon momentum reaches the charge Fermi momentum (of the 
attractive model) = kp = A;|.(att). We shall find power law behaviour for 
charge density correlation functions at momentum transfer q 2kF- 

Finally, the phase shift is 



+ / dk'aAK' 

J\M>B„ 



A,hl AfeN _ 1 1 „ j ^ V ^ \U\ )^ \2^^ \U\ 



^ "|r(i + zA^)r(|-zA^ 

1, fr(i-^ ^)r(| + zA^ 

Aft-A'\ -p / 1 •A'»-A 




L (.V l«l / V2 \u\ 

with the first term corresponding to the triplet phase shift in eq(58). 
Spin excitations 

The charge excitations we studied were gapless at half filling and away 
from it since they essentially involved moving pairs around. The construction 
of spin excitations, on the other hand, requires breaking of pairs and therefore 
an energy gap opens up. 

The Triplet. The k 2-strings describe spin singlets, hence the only way to 
generate a spin excitation from the ground state is to break one of the pairs 
creating momenta k2 the real line with the corresponding the quantum 
numbers and ^2. Now M[ = N/2 — 1, hence the state is a spin triplet, 
consisting of two objects, the spin-1/2 spinous of the attractive model coupled 
symmetrically to form a spin-1 state. 

No hole opens in the {/^} set since one slot less is available to the re- 
duced number of pairs. The quantum numbers will still be distributed sym- 
metrically around the origin. As a consequence the only contribution to the 
excitation momentum comes for the real A;, 

AP = ^{m+n2). (150) 

Since there are no holes in the A-distribution, «t(A) satisfies the integral 
equation 

Ca{A) = cf)o{A) - yi^i(A - sinA;i) - yi^i(A - sinAjs) (151) 

Ij Ij 
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As usual, we introduce an excitation density, (Ti(A), via «t(A) = «To(A) + 
■^«Ti(A) and we have 

^i(A) = ^t^(A) + ^^(A), (152) 
where a^^ is the solution of the equation 

jCa^' = -Ki{A - sin (153) 

or, equivalently 

J^a^^ = a^^'^. (154) 

The excitation energy is 

AE'"^h,k2) = Uh) + e.,ik2) , (155) 

with 

e,,(h) = -2tcoskj - At J^" dAa^' (A) i^Re^l - (A - i^-^Y + ^| ~ 
= + / , t?AV?(A){C(A) + ^ + 



'|A| 

and the excitation momentum 

AP'"^{kr,k2)=Ps{kr)+p,{k2) , (156) 

with 

Psikj) = PTiki) - / iAVj^(A') (- - 2arctan ( eFt(^'-^-*'^)') 1 . (157) 

J\A\>Bo 12 V /J 

Finally, the phase shift is 



S'"P = f " dAar(A) \ 0(2(sinA;i - A)) - 2 ^7 ^' \ (158) 
j-Bo [ an \ 



c?TC v/|A|>Bo I V J 2 (in 
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At half filling we have 

1 , /27r 



hence 



a'l(A) = sech f — (A - sin k)] , (159) 

w,x \V\ , /"°° , Ji(») cos(»sin A;) 

= h^-2tco,k + U dp^ WT^- (1^0) 

2 ^0 ]_ _|_ gVlpl 

Jo{p) sm{p sin k) 
1 + e¥W 

From (160) and (161) we conclude that the spinons underlying the spin 
triplet correspond to the charged holons in the repulsive Hubbard model. 
This is further borne out by the phase shift, 



pTik) = k + 2l dp "^^^ (161) 



/"°° 1 / 27r \ 
-27r / JA— sech —(A - sin k) 0(2(sin k2 - A)), 

J-oo \u\ \\'^\ } 

- log <^ -) (—7^ 

i r (1 + r -i^ 



- log <^ ^ f . (162) 

i I r ( 1 -I- smfci -smfc2 \ -p ( 1 _ - smfci -smfc2 \ 

\u\ J \2 \u\ ) ) 

Away from half filling the symmetry is broken but the identification still 
survives. 

The Singlet. We wish to break a pair without changing the spin, hence 
need to introduce an additional /„. We have M[ = N/2 — 1, Mi = 1. As 
before, we have two real k and quantum numbers TCi, n2. It can be seen 
from (198) that = 0. However, the presence of a new A implies that 
the quantum numbers have to be shifted by | with respect to those in the 
triplet, as can be seen from (197). As a result, an extra tt appears in the 
phase shift. Also, we determine the position of the added spin momentum 
to be A = |(sin ki + sin k2). 

We find that the densities «t(A) and cri(A) are the same as in the triplet 
case. So are the excitation energy and momentum. The only difference with 
respect the triplet case appears in the phase shift, 

gsing ^ gtrip _ q^^-^ _ gin ^=2) - TT (163) 

At half filling, we have 

i^(sin k, - sin k2) r (1 - i -^"^'|-r^O r (H 




°° i I 1 + Z|^(sin fei — sinA;2) T (^1 + ^smfc, -sinfca ^ ^ {2 ~ 
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a result we met as the holon-holon phase shift in the repulsive case, We 
must conclude that the spin excitations here are made of the same objects 
that made the charge excitations in the u > 0. 

Combining the singlet and triplet reults we find that the spin S-matrix 
at half filling has the familiar form for spinon scattering, 



^ T[l-^^-^^^^^^}T[l + ^^^^ ( (sin fe, - sin fe,)/^^ + z^pi^ 

^ sinfci-sinfca ^ r (I - i ^i"fcl-p"fc2 ) | (sin ^ - SIVL Ajj) + 



(164) 



Charge-spin excitations 

We consider now excitations where the number of electrons changes. 
Again we add an extra term to the hamiltonian to make apparent the particle- 
hole symmetry. 

To remove an electron we must break a pair. Once the electron is removed, 
we have an unpaired real k left. We have N — 1 electrons and M' = N/2 — 1, 
hence the number of slots available does not change. As there is one less k 
2-string a hole in the {A} sequence will appear. The state is labeled by the 
holon parameter k'^ and the spinon parameter A'*, with the corresponding 
quantum numbers n'^ and J'*, and momentum momentum excitation AP = 

The integral equations for «t(A) and cr^(A) are straightforward leading to 

AE''\k\ A'^) = e,(A;'^) + e,(A'^), AP''\k\ A'^) = p,(k'') + p,(A'^) ,(165) 

Again we find theat an elctron is comopsed of a spinon excitation carrying 
the spin content and a holon carrying the charge. An electron with definite 
momentum p removed from the system will be a superposition of the type 
of state just constructed subject to AP*''*(A;'*, A'*) = p, the spread being a 
measure of its lifetime. 

The spinon-holon phase shift is 



S''^{k^,A^) = Q{2{smk^ - A^)) 

dK[al{K- K^) + af{K)] |0(2(sin k^ - A)) - 2^^4^ 
Bn dn 



i: 



e - 2^ (166) 
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+ / JA[<(A;A'^) + (7f (A)] i2arctan feF 

^|A|>Bo I V J 2 an 



|(Aft-sinfc'')^ dp. 



At half filling we have, 
1. From (139) 



2. From (159) 

1 , /27r 



(7^'°°(A;A'^) = ^i? ( ^(A- A'^) I , (167) 
\u\ \ \u\ I 



(rf '°°(A) = --sech f — (A - sin k'')] , (168) 
u \ u ) 

3. From (140) 



•p cosh ( 

4. From (160) 

\V\ ,1, Z""" , Ji(») cos(»sin A;'*) 

eT(^') = V-2^^°«^ +4^ '^^^^^^^ ^ (170) 

z ^0 ]_ _|_ e 2 IpI 

5. From (137) 

■^0 pcosn^-'j-pj 

6. From (161) 

= +2 r,i /°w-"'^j;"^''' . (172) 

7. From (166) 

e = ^ + 2 arctan (^eFT(^''--^'')) (173) 
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Once more we see that an energy gap is present at half filling. 



When we add an electron we introduce a real k describing the antiholon 
and, from (198), it is clear that the number of available slots for the 
quantum numbers increases by one, while the same amount of k 2-strings 
are present. As a result, a hole appears in {A}. The state is labeles by the 
spinon parameter A'* and antiholon parameter k'^. 

The integral equations are straightforward and similar to the previous 
case. However, the value of the integration limit is different, since the number 
of k 2-strings does not change. We have J^g dA(T(A) = J^g^ dA(To(A) = 
N/2L^ and hence a shift in the integration limit 

, dAa'JA) , , 

LiB - Bo) = -'-f,^^ (174) 



2 dS, 



This change affects the values of the excitation energy and the phase shift, 
but not excitation momentum, 

= e,(A;'^) + e,(A) + M + 2// , 

AP'''\k^,A^) = p,{k^) + p,{Ah) , 
S'''^(k^,A^) = @(2(smk^ - Ah)) 



+ dA[al(A;A^) + af(A)] |0(2(sinA;'^ - A)) - 2*^^'^^ 
S-^^ + [ dA[al(A; A'') + af(A)] \ 2 arctan (e 

J\M>Bo \ 



dn 

lf|(Aft-sinfcS) 



Turning to half filling, the results (167-173) apply here too. In particular, 
we find that the phase shift is the same 

cs,h A"^'^ 

"n=l — "N=L 

reflecting the appearance of the charge SU{2)^ as the relation between the 
energies 
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reflects the particle hole symmetry. 



We have completed a detailed discussion of the elementary excitations 
and their varoius characteristics. These are the low lying solutions of the 
Bethe-Ansatz equations. To compute the free energy, however, we need the 
complete set of solutions as given by the string hypothesis. 

The Thermodynamics of the Hubbard Model 

The thermodynamics is derived [5] by the method [22] we already em- 
ployed in Lecture 3. We shall merely outline the main steps and write down 
the answer. 

Consider the thermodynamic potential, = E — AN — TS , with S the 
entropy and A the (external) chemical potential. One calculates it summing 
over all energy eigenstates, or equivalently integrating over all allowed solu- 
tion densities. Denote by p{k),an{K) and cr^(A) the distribution functions 
of A;, A("\ A'(") respectively, where A^") is the real part of the A ra-string and 
A'(") is the real part of the k — A ra-string, and by p^{k)^ a^{A) and cr^'*(A) 
the corresponding hole-distributions, see Appendix. Further define: 

C{k) = p\k)lp{k), 7/4A) = ^^(A)/^4A), 7/;(A) = <'^(A)/<(A). 

The same steps that led to the thermodynamic equations in Lecture 3 
lead now to the following expression for the thermodynamoc potential, 

/TT /"OO 
Po{k)\n{l + C{k))dk + / (7„(A)ln(l -7/i(A))JA 
-TT J — OO 

with Eo,(To{A), po{k) being the ground state energy and densities at half 
filling, and the functions Ci^) ^^'^ V'i-i^) ^^^^ determined from the following 
set of coupled integral equations. 



ln7/„ = G'[ ln(l + 7/„+i) + ln(l + 7/„_i) ] tc = 2,3,... 
luT/; = G'[ln(l + 7/;+,)+ln(l + 7/;_,)]r^r^ = 2,3,... 

luT/i = G'[ ln(l + 7/2) - / ^(A - sinA;)ln(l + C"^)cosA; JA;] 

J —TT 

In T][ = G[ ln(l + rj'^) - T ^(A - sin k) ln(l + () cos k dk] 

J —TT 

2 COS k 1 I'K r 4 w 

ln( = — \ — / JA sech — (A — sin A;) — — (sgn [/) i?e^/l — (A — yz)^ + In 

± U J — QO U 



r ^ ° ' V ^ 4 ' 1 + 7/1. 
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where G is the integral operator, Gf{K) = ^ / ^^^ij^ZoZITT /(A')- The asymp- 
totic conditions are, 

\a.rin — > — )• oo (177) 
InT/j^ — )• , as — )• oo 



We leave it as an exercise to show that in the limit T — )• and no magnetic 
field the equations collapse, depending on the the sign of [/, to their respective 
ground state equations: eqs(13,14) in the case of repulsion, eq(lll) in the 
case of attraction. If the zero temperature limit is taken in the presence of a 
magnetic field one obtains the corresponding magnetization equations. 

In the repulsive case, 

1 rB(h) 
p[k) = \-cosk dAa{A)Ki{smk — A) 

27r J-B(h) 

aiA) = / dkpo(k)KUsmk- A)- JAV(A')i^2(A - A'), 

J-Q J-B(h) 

In the attractive case we have 



1 1 rQ(f^) rB 

aiA) = -Re { \ - / KUA-smk')p(k')dk'- JAV(A')i^2(A - A') 

1 fB 

p[k) = \-cosk dAa{A)Ki{smk — A) 

2'ir J-B 

The integration limits are now determined by the imposed magnetic field 
h. In the repulsive case it mainly affects B = B{h) < oo ( is mainly 
determined by the density though there is also a weak dependence on h.) The 
situation is reversed in the attractive model where B is mainly determined 
by the density and Q mainly by the magnetic field. 

Unlike the magnetization equation discussed in lecture 3 the limits are 
imposed symmetrically and therefore analytical results are available only in 
some limits. For extensive numerical work see [17] and [4] in the repulsive 
case, and [18] [19] in the attractive case. 
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The thermodyanmic equations can now be analysed in a manner discussed 
in Lecture 3 to determine the behavior of the model in the infra red. The 
physics is determined by the gapless excitations: charge excitaions in the 
attractive model, spin excitations in the repulsive model at half filling and 
both charge and spin away from it, and will flow to a fixed point accordingly. 

Alternatively, the nature of the fixed point can be identified using Bethe- 
Ansatz finite size calculations combined with methods of conformal field the- 
ory [20] [21]. But that would be the subject of Lecture 5. 

Appendix 

Here we discuss the different species of {%} and {A^} that appear in 
the attractive and repulsive cases. We can classify the eigenvalues in the 
following categories: 

1. Real kj. They have associated quantum numbers rij. 

2. k 2-strings. We may have solutions with pairs of complex A;, ka = 

To be solutions of the Bethe-Ansatz equations, eqs(l,2), in the ther- 
modynamic limit we must have [6], 

-n (4") ± ^x^^^) = A. T ^ J + O (e-"''^) (179) 
This is valid only if we can neglect O ^e"''" terms, where 

7/i±) = xt^ ± j^Im I 2arctan (| (sin (4±) ± ixt^) - A^j] Uo) 

Therefore, we must have xit"^ ^ 0. Let's prove this point. That is, the 
values described by the string hypothesis (180) are indeed solutions of 
the eigenvalue problem. We will start with the spin equation (2). Let's 
write 

sin k± = sin (k it ixa) = A + e=pz— , e^A (181) 

where A belongs to the {A^} set and, for simplicity, we let all the A's 
to be real. Then (2) becomes 
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which is again spin equation, but for N — 2 particles see eq(79). 
We now turn to the charge equation, (1). After substitution, we get 



e — I- , 



T-r Jl-v ~ Sill yri'a. "T t-Xa) ~ , . 

n T^^ , , (183) 



e 



A^ 


- sin 


+ iXa) 




A^ 


- sin 


+ ^Xa) 




A^ 


- sin 


- ^Xa) 




A^ 


- sin 







n 7 : 7 r , » 



Dividing (184) by (183), and taking logarithms we get 
XaL ~ - log e + log - + ^ /m (-2 arctan (sin + ixa) - A^))(185) 

Notice that log ^ <^ L. Therefore, in order for (185) to be satisfied, we 
must have 

e ~ O (e-''^) (186) 

with r/a given by (180). 

When we neglect the exponential term, (179) has two possible solutions 
(a) 

U 



± ixlt^ = arcsin (a„ T i j) , < < ^ 



cos kI^^ 
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This gives 

sin (^kI^^ ± ixlt^) = sin k^^^ cosh xlt^ ± i cos k^^^ sinh xlt^ = A„ =F 

(188) 

Since cos K^^^ > 0, 

xt^>o^u<o 

(189) 

(b) 

± ixi^^ = ^ - arcsin (a. t i j) , ^ < < y ^ cos <(±) < (IL90) 
This gives 

sin (^k'J^^^ ± ixlt^) = sin k'J^^^ cosh ± z cos k^^^^^ sinh xlt^ = An T i j 

(191) 

Since cos < 0, 

> > (192) 

We conclude that the forms of the k 2-string that is a solution of the 
problem depends on the sign of the interaction. 

Therefore, the k 2-strings are of the form 



i) u> 0. 



ii) u < 0. 



= TT — arcsin(A' — i—) 
u 

k+ = TV - arcsin(A' + i-) , (193) 



K = arcsin(A' - i^) 

= arcsin(A' + z-^) . (194) 
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Notice that each pair is described by only one parameter, A^. The 
quantum numbers associated to the pairs will be denoted 

3. Quartet, k — K 2-string. One of the possible excitations in the attractive 
case involves a group of four complex ka such that 



kl = arcsin(A" + z y) , 

kl = TT - arcsin(A") , 

kl = arcsin(A") , 

kl = arcsin(A:: - ) (195) 



parametrized by a single variable A" and the quantum number associ- 
ated is 

4. Real A^. The corresponding quantum numbers are denoted /„, as in 
the ground state. 

5. A 2-string. It will only appear in the elementary excitations for u > 0. 
We have 

A«=A«±i^. (196) 

The pair is parametrized by the real and has quantum number 
associated to it. 

In the thermodynamic limit, where the string hypothesis is valid, the 
objects enumerated are allowed solutions of the Bethe-Ansatz equations, and 
the sets {A^}, {A"}, and {A^} are subsets of {A^}. 

We wish to count the number of slots available for each configuration. 
Denoting, 



M 


= number of real A^, 


M 


= number of A^, A 2-strings 


M' 


= number of A^, k 2-strings 


M" 


= number of A", quartets 


M'" 


= M' + 2M", 


N - 2M'" 


= number of real kj, 


Ni 


= M + 2M + M'", 




= N-N^. 
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we have (assuming L of the form, X = 4i/ + 2, integer 

_ J integer (M + M + M' + M" even), \nj\<Ll2 
"■^■"Ih.o.i (M + M + M' + M" odd), \n-\<{L-l)l2 ^ ' 



< ^(iV-2M"'-(M + 2M)-l) , /„ 
\J^\ <'^{N - 2M'" - 2M - M - I) , J« 
\r^\<-{L-N + 2M"'-M' + 2M"-l) , /„ 



|/''| < -(X - iV + 2M'" - 2M' + M - 1) , / 



integ 


;er 


{N- 


■ M 


odd) 


h.o.i 




{N- 


■ M 


even) 


inte^ 


fer 


[N- 


- M 


odd) 


h.o.i 




{N- 


- M 


even) 


integ 


;er 


(L- 


(N 


-M') 


h.o.i 




(L- 


(N 


-M') 


integ 


;er 


(L- 


{N 


-M") 


h.o.i 




(L- 


(N 


-M") 
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